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Polytechnical University, Sankt-Petersburg 195251, Russia 


Abstract 

The exact nonperturbative confining solutions of the SU(3)-Yang-Mills equations 
recently obtained by author in Minkowski spacetime with the help of the black hole 
theory techniques are analysed and on the basis of them the gluon propagator corre¬ 
sponding to linear confinement at large distances (small momenta) is constructed in 
a nonperturbative way. At small distances (large momenta) the resulting propagator 
passes on to the standard (nonperturbative) gluon propagator used in the perturbative 
quantum chromodynamics (QCD). The results suggest some scenario of linear con¬ 
finement for mesons and quarkonia which is also outlined. As a consequence there 
arises a motivation for studying the relativistic bound states in the above confining 
SU(jV)-Y ang-Mills fields. This possiblity is realized for N = 2, 3,4 with the aid of 
the black hole theory results about spinor fields on black holes with a subsequent ap¬ 
plication to the charmonium spectrum in the most important physical case N = 3. 
Incidentally uniqueness of the confining solutions is discussed and a comparison with 
the nonrelativistic potential approach is given. 
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1 Introduction and Preliminary Remarks 

As soon as quantum chromodynamics (QCD) was proposed as the main candidate for the 
theory of strong interactions (I] at once there arised the question about the confinement 
of quarks within the framework of QCD. In the late seventies of XX century the main ap¬ 
proaches to solve the problem were formed (see, e.g., review of Ref. 0 ) and they actually 
remain the same ones up to now. For this purpose miscellaneous techniques were elabo¬ 
rated, for example, strong coupling expansions, lattice approach, instanton improvement of 
perturbation theory, nonrelativistic potential approach and so on. It should be noted, how¬ 
ever, none of the mentioned various directions has so far led to a generally accepted theory 
of quark confinement. 

In this paper we would like from another side to have analysed one of the possible ap¬ 
proaches. The question is about a nonperturbative modification of gluon propagator which 
might correspond to linear confinement between quarks at large distances. The very sim¬ 
ple idea of modifying the mentioned propagator arises when considering the naive Fourier 
transform for the power potentials of form r x (for more details, see e.g. Ref. 0 ). Then 
at A = 1 (linear confinement) the conforming Fourier transform (propagator) is of order 
|k| -4 in momentum space, while the case A = —1 (Coulomb potential) gives the standard 
gluon propagator ~ k 2 . All the attempts to obtain the necessary behaviour, however, 
for example, by summing a selective infinite set of perturbation diagrams with using the 
Dyson-Schwinger equations for the propagator failed 0. It is clear why: it is impossible 
to get anything nonperturbative like confinement by perturbative techniques. Some new 
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possibilities in this direction were connected with lattice theories (for more details, see, 
e.g., Ref. |j3) and references therein) but the results here are mainly of qualitative character. 

To our mind, from the very outset the problem should be considered on the basis of 
the exact nonperturbative solutions of the SU(3)-Yang-Mills equations modelling quark 
confinement which, in what follows, we shall call the confining solutions. Such solutions 
will be supposed to contain only the components of the SU(3)-field which are Coulomb-like 
or linear in r, the distance between quarks. In Ref. j4j a number of such solutions have been 
obtained and the corresponding spectrum of Dirac equation describing the relativistic bound 
states in those confining SU(3)-Yang-Mills fields has been analysed. Further in Refs. Eli 
M\ the results obtained were successfully applied to the description of the quarkonia spectra 
(charmonium and bottomonium). In its turn, the mentioned description suggests that linear 
confinement is (classically) governed by the magnetic colour field linear in r and, as was 
mentioned in Refs. El EH, one can try to modify the gluon propagator nonperturbatively 
at quantum level for to generate the mentioned magnetic colour field at classical level. One 
part of the present paper is just devoted to it. It should be noted, however, that all the 
main features of such a modification may occur already within quantum electrodynamics 
(QED) that should be not surprising because, as is historically known (see, e.g., Ref. □), 
the standard gluon propagator is in fact the slightly modified photon one of QED. Under 
the circumstances we shall conduct our considerations mainly within QED but incidentally 
making remarks to generalize the results obtained to the QCD case. Considerations of the 
first part of paper (Sections 1-4) inevitably lead to the task of a more thorough analysis of 
the Yang-Mills and Dirac equations derived from QCD lagrangian which is realized in the 
rest of paper. 

Let us introduce some notations. Further we shall deal with the metric of the flat 
Minkowski spacetime M that we write down [using the ordinary set of local rectangular 
(Cartesian) (x, y, z) or spherical (r, 0, <p) coordinates for spatial part] in the forms 

ds 2 = g^dx^ ®dx v = dt 2 — dx 2 — dy 2 — dz 2 = dt 2 — dr 2 — r 2 (dd 2 + sin 2 ddp 2 ), (1.1) 

so the components <y /( „ take different values depending on the choice of coordinates. Be¬ 
sides we have S = | dct(g lw )\ = (r 2 sin$) 2 in spherical coordinates and the exterior dif¬ 
ferential d = dtdt + d x dx + dydy + d z dz or d = dfdt + d r dr + d$dd + d v dip in the corre¬ 
sponding coordinates. We denote 3-dimensional vectors by bold font so X = (f, x, y, z) = 

(t, r), k = (k 0l ki,k 2 ,k 3 ) = (k 0 ,k) with r 2 = x 2 + y 2 + z 2 , |kj = yjk\ + k^ + k%. 

The Fourier transform T>(A:) of some function 4>(X) is formally defined by the relations 

(kX = kot — k\x — k 2 y — k% z) 


$(fc) = f exp(ikX)<S>(X)d 4 X = F[$] , 


Jm 



( 1 . 2 ) 


but it is treated in the sense of the theory of generalized functions (distributions) (see, 
e.g., Refs. 0) and we denote d 4 k = dkodk\dk 2 dk^, d 4 X = dtdV, dV = dxdydz or 
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dV = '/SdrdOdif while for the generalized 5-functions we use the notations 5(X) = 
6(t)6(x)S(y)d(z), 5(r) = 6(x)S(y)6(z). Other mathematical results necessary for our con¬ 
siderations are gathered in Appendices A, B,C, D, E. 

Throughout the paper we employ the system of units with h = c = 1, unless explicitly 
stated otherwise. Finally, we shall denote L‘A, F) the set of the modulo square integrable 
complex functions on any manifold F furnished with an integration measure while Lf) (F) 
will be the n-fold direct product of I^AF) endowed with the obvious scalar product. 


2 Confining Solutions of the Maxwell and SU(3)-Yang-Mills 
Equations 

2.1 Black Hole Theory Techniques 

For obtaining a set of the confining solutions within the given paper we shall employ the 
techniques used in Refs. 0 for finding the U(A r )-monopolc solutions in black hole physics 
and the essence of those techniques consists in systematic usage of the Hodge star operator 
(see Appendix A) conforming to metric (1.1). As is known, such a metric can be obtained 
from the Schwarzschild black hole metric when the black hole mass is equal to 0. 

Really, if writing down the Yang-Mills equations (B. 3) in components then we shall 
be drowned in a sea of indices which will strongly hamper searching for one or another 
ansatz and make it practically immense. Using the Hodge star operator as well as the rules 
of external calculus makes the problem quite foreseeable and quickly leads to the aim. 

As was remarked in Appendix B, the sought solutions are usually believed to obey an 
additional condition and as the latter one in the present paper we take the Lorentz condition 
that can be written in the form 

div(A) = 0 , (2.1) 

where the divergence of the Lie algebra valued 1-form A = A^T a dx^ is defined by the 
relation 

div(A) = ±=d„(V6g^A u ) . (2.2) 


2.2 Electrodynamics 

We proceed from the second pair of Maxwell equations {B. 5) 


d * F = J (2-3) 

with F = dA, A = A^dxJ' and the Hodge star operator * is defined, for example, on 2- 
forms F = F ltl/ dxJ l A dx v in Minkowski spacetime M provided with a pseudoriemannian 
metric g jW (1.1) by the relation (see Appendix A) 

F A *F = - g^g ua )F^F^V6 dx 1 A dx 2 • • • A dx A 


(2.4) 
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in local coordinates x = (x 11 ) while J = jy* (dx 11 ) with a 4-dimensional electromagnetic 
density current j = jyd/xJ 1 . Let J = 0 and we shall search for the solution of (2.3) in the 
form A = A t (r)dt + A,p(r)dip. It is then easy to check that F = dA = —d r A t dt A dr + 
d r A v dr A dip and since *(dt A dr) = — r 2 sin tid'd A dip, *(dr A dip) = — s\t$dt A dr? we 
get = ?’ 2 sin A t dd A dip — ^Pdd r A v dt A cW. From here it follows that Eq. (2.3) 
yields 

dr (r 2 dr At) = 0, d 2 r A v = 0, (2.5) 

and we write down the solutions of (2.5) as 

A t = —b A , Atn = br + B (2-6) 

r 

with some constants a, b, A, B parametrizing solutions (further for the sake of simplicity let 
us put a = 1, b = 1 GeV, A = B = 0). 

To interpret solutions (2.6) in the more habitual physical terms let us pass on to Carte¬ 
sian coordinates employing the relations 


ip = arctan(y/x), dip 


dl P , , dip 

—dx + —dy 
ox ay 


(2.7) 


which entails 

ry rx 

Apdip = - 5 ——^dx + —r——^dy ( 2 . 8 ) 

x z + y z ar + 

and we conclude that the solutions of (2.6) describe the combination of the electric Coulomb 
field with potential <f> = 1/r and the constant magnetic field with vector-potential 


A — (A x ,Ay,A z ) — I 


ry 


rx 


x A + y 2 ’ x 2 + y 


:,2 ’ 


0 = - 


sm ip cos ip 
sin d ’ sin id 


,0 


(2.9) 


which is linear in r in spherical coordinates and the 3-dimensional divergence divA = 0, 
as can be checked directly. Then Eq. (2.3) in Cartesian coordinates takes the form 


A4> = 0, rotrotA = AA = 0 (2.10) 

with the Laplace operator A = d 2 + d 2 + d 2 . At last it is easy to check that the solution 
under consideration satisfies the Lorentz condition (2.2) 

div(A) = -^j=d tl (y/8g ,xv A v ) = 0 . (2.11) 


2.3 SU(3)-'Yang-Mills Theory 

Now Eq. (2.3) should be replaced by the Yang-Mills equations (B. 3) 

d * F = g(*F A A — A A *F) + J (2.12) 

for SU(3)-field A = A^dx^, Ay = where the matrices T a form a basis of the Lie 

algebra of SU(3) in 3-dimensional space, a = 1,..., 8 and further let us put T a = \ a , where 
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X a are the Gell-Mann matrices (see Appendix B). After this we search (at J = 0) for the 
solution of (2.12) in the form A = A t (r)dt + A^(r)d<£> with At >tp = Aj^Ay + Af^As. 
Evaluating dF = <7.4 + gA A A it is easy to check that the right-hand side of (2.12) is 
equal to zero and to gain the sought solution in the form (which reflects the fact that for any 
matrix T from SU(3)-Lie algebra we have Tr T = 0) 

A t + ^ A t = -^ + At,-A| + i=A t 8 = -^ + A 2 ,--^=A 8 = S±±*L- {Al + A 2 ), 

= bir+Bi, —A^ + —= 6 2 r+.B 2 , ——= — (6i+6 2 )r — (B\ +i? 2 ), 

(2.13) 

where real constants aj , Aj , bj , Bj parametrize the solution, and we wrote down the solu¬ 
tion in the combinations that are just needed to insert into the corresponding Dirac equation 
(see Section 5). From here it follows 

At = [(®2 — Qj\)/r + A\ — A 2 ]/2, A® = [Ai + A 2 — (ai + o 2 )/r]\/3/2 , 

A^ = [(&! - b 2 )r + B 1 - B 2 \/ 2, A« = [(h + b 2 )r + B X + B 2 }V 3/2 (2.14) 

and practically the same considerations as the above ones in electrodynamics show that the 
given solution describes the configuration of the electric Coulomb-like colour field (com¬ 
ponents A t ) with potentials <b 3 . <h s and the constant magnetic colour field (components 
A<p) with vector-potentials A 3 , A 8 which are linear in r in spherical coordinates with 3- 
dimensional divergences divA 3 =divA 8 = 0 while the Eq. (2.12) is easily transformed 
into the Eqs. (2.10) with an obvious modification. It is also simple to check that the solu¬ 
tion under consideration satisfies the Lorentz condition (2.2). 


3 Linear Confinement in QED 

Now we can investigate how the photon propagator should be modified if in the real world 
the interaction between two charged particles would not be classically described only by 
the Coulomb law but it would also include a constant magnetic field linear in r, the distance 
between particles, so the given field would obey the Maxwell equations. Let us briefly 
recall the scheme in accordance with that the standard photon propagator is obtained (see, 
e.g., Ref. Q). 

3.1 Standard Photon Propagator 

In Cartesian coordinates the Eq. (2.3) for A = A^dx^ at J = 0 with Lorentz condition 
div(A) = 0 reduces to the system 

□a m = o 


(3.1) 
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with the d’Alembert operator □ = d 2 — A. Then one constructs the fundamental solution (a 
Green function) of the system (3.1) (the photon propagator at classical level) as the matrix 
with elements D flI/ (X) obeying the system 

ga^DfiuiX) = -g av 5(X) . (3.2) 


Further one uses the fundamental solution of the d’Alembert operator [i.e., the solution of 
the equation OK = 5(A)] in the form 


K(X) 


1 

47r 2 i(A 2 - iO) 


(3.3) 


with quadratic form X 2 = t 2 — r 2 (the exact definitions concerning the generalized func¬ 
tions connected with quadratic forms can be found in Refs. @), so that 


D IW (X) 


-g^K(X) = 


_ 9fu' _ 

47r 2 i(A 2 - iO) ' 


(3.4) 


After this the photon propagator at quantum level is obtained as the Fourier transform for 
D^ W (X), namely 

D, u (k) = F[D^(X)] = , (3.5) 

since K = F[K] = 1 /(A: 2 + /.()) with quadratic form k 2 = k$ — k 2 — k 2 — k 2 (for more 
details see Refs. @). There arises the question why among a large set of the mathematically 
possible fundamental solutions (Green functions) for Eq. (3.2) one chooses just the one of 
(3.4). The answer can be based only on physical considerations. 

Indeed, let us take a point particle with a charge e moving with a velocity v = v(i). 
Then, as is known (see, e.g. Ref. Col), the 4-dimensional density current of such an object 
is j = j^dx 11 ' = ed(r)(dt + vdr), dr = (fix. dy, dz). Under the circumstances the particle 
will generate the electromagnetic field with potential A = A^dx^ which should be obtained 
by the contraction of j with a fundamental solution D^(X) of Eq. (2.3) or, that is the same, 
of Eq. (3.1), namely 


A„(X) = 4irJ D^(X - X')j v {X')d*X' . 


(3.6) 


When choosing D IW (X) equal to that of (3.4) we obtain, for example, for electric potential 
of the field generated [with replacing t —» it in the integral over f in (3.6)] 

r°° p 

MX) = *(X) = - / = - , (3.7) 

7r J —oo t + r z r 

that is, the Coulomb law. Analogously, other diagonal components of D fiv (X) of (3.4) 
give, for example, at v = const the vector-potential A for a Coulomb-like magnetic field 
generated by the point charged particle when its moving (for more details see Ref. llOl ). 
All of that corresponds to experimental data and, as a result, the choice of D jlv (X) in the 
form (3.4) reflects the real situation in our world. 
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3.2 Electrodynamics with Linear Confinement 


Let us now explore how the photon propagator should be modified if experiment would say 
to us that when its moving a point charged particle generates a constant magnetic field linear 
in r, the distance from the particle, additionally to the mentioned Coulomb-like fields. 

Under the circumstances we should use the property of any fundamental solution that 
the latter is determined only to within adding any solution of the conforming homogeneous 
equation. As we have seen above, components of A from (2.9) are the solutions of the 
Laplace equation. On the other hand, as is known, the Coulomb potential is also the solu¬ 
tion of the Laplace equation at r ^ 0 and besides it is a fundamental solution of the Laplace 
operator |9). Consequently, for to obtain the necessary modification of 3-dimensional pho¬ 
ton propagator we should add components A x or A y of (2.9) to the Coulomb fundamental 
solution for the conforming components of propagator. To pass on to a 4-dimensional prop¬ 
agator let us recall that the Coulomb fundamental solution of the Laplace operator and 
the fundamental solution (3.3) of the d’Alembert operator are connected by the so-called 
method of descent (see, e.g.. Ref. liTTl l which is in essence expressed by the relation 

r°° 1 

4vr / K(X)dt = — , (3.8) 

4-00 r 


so that to obtain the necessary modification of 4-dimensional propagator we should add 
some suitable solutions of the d’Alembert (wave) equation □/(X) = 0 to K{X ) for the 
corresponding components of propagator (3.4). As the latter ones we should take functions 
(tA x )/ (87 t 2 ) or (tA y )/(87 t 2 ) [where the factor 1 /( 87 r 2 ) is introduced for the sake of further 
convenience] with A x or A y of (2.9). The given functions obviously satisfy the d'Alembert 
equation. Then the sought photon propagator will have the same components as in (3.4) 
except for the cases fi = u = x or y where the components will be 


D xx (X) = 


Dyy(X) = 


1 


T 


tA q 


1 


try 


47 r 2 i(X 2 — iO) 87 t 2 47 r 2 i(X 2 — iO ) 8ir 2 {x 2 + y 2 ) 


1 


+ 


t,A 


y _ 


1 


+ 


trx 


4ir 2 i(X 2 — iO) ' 8ir 2 47 r 2 *(X 2 — *0) ' 87 v 2 (x 2 + y 2 ) ^ ^ 

Under this situation when its moving a charged particle might generate an additional con¬ 
stant magnetic field according to the relation (3.6). To specify it let us recall that in elec¬ 
trodynamics in any constant magnetic field is connected with a finite motion of charged 
particles. Let us suppose, for instance, that the particle accomplishes a finite motion within 
a finite (though perhaps large enough) time in such a way that the velocity projections 


J x,y 


(t) are some odd functions of time. Then one may consider that 


jit- t')v x ^ y it')dt' ~ J t' v Xt y{t')dt' ~ C = const , 


(3.10) 


and according to (3.6) there appears some constant magnetic field ~ A of (2.9) linear in 
?% the distance from particle. So indeed under the certain conditions we could observe the 
mentioned magnetic field corresponding to the propagator described above. 
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3.3 Momentum Representation 

To get the necessary propagator at quantum level we should carry out the Fourier transform 
of the just found propagator. So long as for any natural m we have (see Refs. f9l l 

F[t m ] =2Tr(-i) m 6 im \k 0 ) (3.11) 

with m-th derivative of h-function, then really everything reduces to the Fourier transforms 
for functions A x ,A y of (2.9). Since the latter ones are only locally integrable the Fourier 
transforms should be understood in the sense of the theory of generalized functions 191 . 
namely, through analytical continuation of suitable integrals. Let us find, e.g., F\A X \. In 
accordance with (1.2) we shall have 

F[A X ] = [ A x exp(-ikr)dV = 


10 


) rn rZn 

r 2 dr / exp(— ik^r cos d)dd / exp [—ir sin d(k\ cos ip + ^2 sin ip)\ sin pdip . 

Jo Jo 

(3.12) 


Using the relation of Ref. d 

r 2n 

exp(a cos x + b sin x) 


>o 


sin x 
cos x 


dx = 


IT 


\J a 2 + b 2 


Ii(\J a 2 + b 2 ) 


2b 
2 a 


(3.13) 


with the modified Bessel function I\{z) = where Ji(z) = —J\{—z) is the stan¬ 

dard Bessel function, we can rewrite (3.12) as 

F[A X \ = ~ j7n ^ 2 _ f r 2 ^ r / exp(— ik^r cos d) J\ (r sin dJk 2 + k 2 )dd . (3.14) 

\]k 2 + k 2j o J o 

Further replacing cos d = x and employing the formula of Ref. ED 

/“ COS ( ^ q2 ~f ) - J l {cx)dx = — [cos (ab) - cos(aV / 62~+c 2 )] , (3.15) 

Jo da 2 - x 2 ac 


we get 


F[A X \ = ^2 + 1.2 J Q rlcosihr) - cos(r|k|)]dr . 


At last, using the relation with the Euler T-function of Ref. D3 

roo 

/ x 9_1 cos {mx)dx = 

Jo 


r (g) nq 

-cos — , 

mi 2 


(3.16) 


(3.17) 


holding true at 0 < q < l,m > 0, we analytically continue the right-hand side of (3.17) 
over all admissible values q, m which permits [at q = 2 with F (2) = 1] to write down 


F[A X ] = - 


47riAi2 
k 2 + 


1 

kl 


1 


47rifc2 
|k| 2 k\ ' 


(3.18) 
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Analogous consideration yields 


F[Ay\ 


Amk\ 
|k| 2 k\ ' 


(3.19) 


After this, employing the relations (3.9) and (3.11) we finally come to the conclusion that 
the sought modification of photon propagator in momentum space will have the same com¬ 
ponents as in (3.5) except for the cases // = v = 1 or 2 where the components will be 


Dxx{F) — 


Fyy(k) 


1 


k 2 + iO 


1 


+ 


k 2 S\k 0 ) 

| k | 2 fcg ’ 

ki5'(k 0 ) 


k| 2 fc| 


wv "' k 2 + iO 

and the generalized function 5' (ko) acts according to the rule 

/ ${ko)S'(k 0 )dk 0 = — 3>'(0) . 


(3.20) 


4 Linear Confinement in QCD 

4.1 Motivation 

It is clear that considerations of the previous section are not confirmed experimentally in 
electrodynamics - there exist no elementary charged particles generating a constant mag¬ 
netic field linear in r, the distance from particle, the given field obeying the Maxwell equa¬ 
tions. 

Another matter is quantum chromodynamics. The analogue of charge here is colour. 
The group U(l) and the Maxwell equations are replaced by SU(3) and the Yang-Mills equa¬ 
tions but, as we have seen in Section 2, both the Maxwell equations and the Yang-Mills ones 
possess the confining solutions. Though quarks can unlikely be considered classical parti¬ 
cles, after all, they accomplish a finite motion within a region with character size of order 1 
fm = 10 ” 13 cm and, as the explicit form of modulo square integrable solutions of the Dirac 
equation (5.4) in the confining SU(3)-field (2. 13)— (2. 14) (relativistic bound states) shows 
(see Refs. HI El [6l|35) and Section 7), the j-th colour component for the system of two 
quarks (e.g., for quarkonia) T'j ~ r aj e^^ jr with ay > 0 , f3j = ^/Xq — lu? + g 2 b 2 > 0 , 
where u = J2 Wj is an energetic level of system, 612 are the parameters of linear interaction 
from the solutions (2. 13)— (2. 14), 63 = — (61 + b 2 ), r is a distance between quarks and is 
proved to be markedly different from zero only at r ~ 1/Pj ~ 0.04 fm (see Refs. @s 6 35i| 
and Section 9), i.e., we deal with linear confinement of colour. Just the magnetic colour 
field defines the latter through the coefficients Bj. As a result, there are certain grounds to 
consider the qualitative physical picture from the previous section to occur just within QCD 
and the gluon propagator should be modified. The necessary modification can be realized in 
the same way as is done when deriving the standard gluon propagator (see, e.g., Ref. 0 ), 
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i.e., through multiplying the propagator (3.5) [where the modification (3.20) is implied] by 
the factor 5 ab with a, b = 1 ,..., 8 


D%[k) = d ab D flv (k) . (4.1) 

Under the cicumstances the gluon propagator obtained will be able to lead to linear 
confinement at large distances (small momenta) while at small distances (large momenta) 
we can omit the additional addenda of order | k | 3 in (3.20) and the resulting propagator will 
pass on to the standard gluon one used in the perturbative QCD. 

It should be noted that during all the considerations we in fact dealt with the so-called 
Feynman gauge (a = 1) but, as is not complicated to see, it is easy to generalize the results 
to an arbitrary a-gaugc since this generalization concerns only the standard part of the 
propagator obtained and has been repeatedly discussed in literature (see, e.g., Ref. [|71 and 
references therein). 

Finally there are four important remarks. 

Firstly, the fact is that the notion of propagator makes no sense for general nonlinear 
equations such as the Yang-Mills ones (2.12). If restricting, however, to the SU(3)-fields 
taking the values in the Cartan subalgebra of the SU(3)-Lee algebra (see Appendix B), 
i.e. in the subalgebra generated by the matrices A 3 and As, then the equations ( 2 . 12 ) (at 
J = 0 ) become linear since the right-hand side of ( 2 . 12 ) is equal identically to zero for 
such field configurations (and for those gauge equivalent to the latter). The confining solu¬ 
tions (2.13)—(2.14) are just of the mentioned class. As a consequence, our modification of 
propagator holds true just in the latter set of SU(3)-Yang-Mills fields but the standard gluon 
propagator is tacitly supposed to correspond to the given class as well because it conforms 
to the Coulomb-like part of the solutions (2.13)—(2.14). 

Secondly, one should say a few words concerning the nonrelativistic confining poten¬ 
tials often used, for example, in quarkonium theory (see, e.g., Ref. Cl)- The confining 
potential between quarks here is usually modelled in the form a/r + br with some constants 
a and b. It is clear, however, that from the QCD point of view the interaction between quarks 
should be described by the whole SU(3)-field A /t = A“T a , genuinely relativistic object, the 
nonrelativistic potential being only some component of Af surviving in the nonrelativistic 
limit when the light velocity c —> 00 . 

As has been mentioned in Refs. EH (see also Section 8), however, the U(l)- or 
SU(3)-field of form .4'/ = Br 7 , where B is a constant, may be solution of the Maxwell and 
Yang-Mills equations (2.3), (2.12) (at J = 0) only at 7 = —1, i. e. in the Coulomb-like 
case. 

As a result, the potentials employed in nonrelativistic approaches do not obey the 
Maxwell or Yang-Mills equations. The latter ones are essentially relativistic and, as we 
have seen, the components linear in r of the whole A fl are different from A t and related 
with (colour) magnetic field vanishing in the nonrelativistic limit. That is why the nonrel¬ 
ativistic confining potentials cannot be assumed as a basis when deriving the modification 
of gluon propagator under consideration. 
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Thirdly, it should be emphasized that the standard gluon propagator (as well as the 
corresponding photon one) is by itself essentially nonperturbative object [a fundamental 
solution (a Green function) of d’Alembert-like system (3.2)] that cannot be calculable by 
any perturbative techniques. Another matter that standard propagator is used in QCD per¬ 
turbation theory for calculating quantum corrections, including for the propagator itself. 

From mathematical point of view we should choose the Green function for the corre¬ 
sponding equations [SU(3)-Yang-Mills ones of (2.12)] which takes into account the nec¬ 
essary boundary conditions - linear confinement which is essentially nonperturbative phe¬ 
nomenon. The standard (nonperturbative) gluon propagator does not satisfy this condition. 
That is why we should modify the mentioned propagator in a nonperturbative way which 
can be done only on the basis of the corresponding exact nonperturbative solutions of the 
SU(3)-Yang-Mills equations. The resulting propagator in (4.1) is essentially nonperturba¬ 
tive one in each summand. 

Fourthly, the structure of the obtained propagator of (3.20) shows that it has rather 
strong infrared singularities at k —» 0. Physical meaning of this is that quarks mainly emit 
and interchange the soft gluons (i.e., those with small /,:), so that gluon concentrations in the 
confining SU(3)-gluonic field are much greater than the estimates given in Section 9 (see 
also Ref. l35l t because the latter are the estimates for maximal possible gluon frequencies, 
i.e. for maximal possible gluon impulses (under the concrete situation of charmonium 
states). It is also clear that just magnetic part of the propagator (3.20) is responsible for 
larger portion of gluon concentrations since it has stronger infrared singularities than the 
electric part. 


4.2 A Scenario for Linear Confinement 

The above results and those of Refs. 0001351 suggest the following mechanism of 
confinement to occur within the framework of QCD (at any rate, for mesons and quarko- 
nia). The gluon exchange between quarks is realized by means of the propagator described 
above. At small distances one may neglect additional addenda of order |k| 3 in the propa¬ 
gator and we obtain the standard gluon one used in the perturbative QCD and, as a result, 
asymptotic freeedom. At large distances the mentioned gluon exchange leads to the con¬ 
fining SU(3)-field of form (2.13)—(2.14) which may be considered classically (the gluon 
concentration becomes huge and gluons form the boson condensate - a classical field) and 
is a nonperturbative solution of the SU(3)-Yang-Mills equations. Under the circumstances 
mesons are the relativistic bound states described by the corresponding wave functions - 
nonperturbative modulo square integrable solutions of the Dirac equation in this confin¬ 
ing SU(3)-field 00 0 (35). For each meson there exists its own set of real constants 
aj,Aj,bj,Bj parametrizing the confining gluon field (2.13)—(2.14) (the mentioned gluon 
condensate) and the corresponding wave functions while the latter ones also depend on //q, 
the reduced mass of the current masses of quarks forming meson 000|3l. It is clear that 
constants aj , Aj , bj, Bj , po should be extracted from experimental data and such a program 
has been just realized in Refs. 00 m for quarkonia (see also Section 9). 
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Finally it should be emphasized that the reason why Nature chose the somewhat differ¬ 
ent propagator for gluons than the one for photons remains obscure within the framework of 
our considerations. To our mind, however, the questions of such a kind may be considered 
only from the cosmological positions, so that, for example, the problems concerning the 
number of quark flavours in QCD or a mechanism of breaking chiral symmetry in QCD can 
be scarcely resolved within the QCD framework in a single-valued way and a cosmological 
approach might be quite useful d. 

5 Confining Solutions of the SU(AO-Yang-Mills Equations 

5.1 Motivation 

The previous section led us to the following problem: how to describe possible relativistic 
bound states in the confining SU(3)-Yang-Mills fields? The sought description should be 
obviously based on the QCD-lagrangian. Let us write down this lagrangian (for one flavour 
and N quark colours) in arbitrary curvilinear (local) coordinates in Minkowski spacetime 


C = m- - g ,l 'g ua )F^Ke, g<v,a<(3 (5.1) 


where, if denoting S(M ) and £, respectively, the standard spinor bundle and ^-dimensional 
vector one [equipped with a SU(A r )-connect ion with the corresponding connection and cur¬ 
vature matrices A = A^dx^ = A a ^T a dx^, F = dA + gA A A = F^ v T a dx^ A dx v , see Ap¬ 
pendix B ] over Minkowski spacetime, we can construct tensorial product S = S(M) 0 £. 
It is clear that T is just a section of the latter bundle, i. e. T can be chosen in the form 
T = (4' |...., 4Ay) with the four-dimensional Dirac spinors ipj representing the j-th colour 
component while 0 = vpt (j° <g> I N ) is the adjont spinor, (f) stands for hermitian conju¬ 
gation, J/v is the unit matrix N x N, <S> means tensorial product of matrices, py is a mass 
parameter, V is the Dirac operator with coefficients in £ (see below). At last, we take the 
condition TV (T a Tp) = KS a b with some real K so that the third addendum in (5.1) has the 
form G(F. F)/ (4 K) (see Appendix A), where coefficient 1 /(4 K) is chosen from physical 
considerations. 

Of course, the most physically important case is that of group SU(3), i.e., three colors 
of quarks, however, it makes sense to have analysed the general SU(A r )-case with arbitrary 
N as well because during a long time there is a firm belief in that considering the limit 
N —> oo to be rather (or even extremely) important for understanding of the real 4D QCD 
with three colors. The spectrum of speculations on this topic ranges from phenomenological 
and lattice approaches (for more details see, e. g.. Refs. ED and references therein) to the 
exotic scenarios based on M -theory and strings (see, e. g., Ref. fT71 and references therein). 

From general considerations the explicit form of the operator V in local coordinates x 1 ’ 
on Minkowski spacetime can be written as follows 



(5.2) 
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where g is a gauge coupling constant, the forms u: a b = uj lia i,dx 11 obey the Cartan structure 
equations de a = u a b A e b with exterior derivative d, while the orthonormal basis e“ = 
in cotangent bundle and dual basis E a = E^d fl in tangent bundle are connected 
by the relations e a (Eb) = 6 b . At last, matrices 7 “ represent the Clifford algebra of the 
corresponding quadratic form in C 4 . Below we shall deal only with 4D lorentzian case 
(quadratic form Q 73 = Xq — x 2 — — xi). For this we take the following choice for 7 “ 


7 


0 






,6 = 1 , 2 , 3 , 


(5.3) 


where 07 denote the ordinary Pauli matrices (see Appendix B). It should be noted that, 
in lorentzian case, Greek indices p,i/,... are raised and lowered with g /w of (1.1) or its 
inverse g 1 ' 1 ' and Latin indices a, b ,... are raised and lowered by = g ab = diag(l,-l,-l,- 
1) except for Latin indices connected with SU(A r )-Lie algebras, so that e“e^ p'"' = r/ ab , 
EaE^g^u = Vab and so on but T a = T a . 

Under the circumstances we can obtain the following equations according to the stan¬ 
dard prescription of Lagrange approach from (5.1) 


, 


(5.4) 


d* F = g(*F A A — A A *F) + gj , (5.5) 

where * means the Hodge star operator conforming to a Minkowski metric, for instance, in 
the form of (1.1), while the source J (a nonabelian SU(A r )-currcnt) is 


J = j“T a * (dx») = *j = *(jZ,T a dx») = *(j a T a ) (5.6) 


where currents 

f = ftdx» = ( 7/i ® I N )T a * dx» , 

so summing over a = 1,..., N 2 — 1 is implied in (5.1) and (5.6). 

When using the relation (see, e. g. Refs. j22j) 

'y c E^LO fia b'y a 'y b = ui fia b'y >1 'y a 'y b = — div( 7 ) with matrix 1 -form 7 = 7 ^dx^ [where div 
is defined by relation (BA)] and also the fact that ( 7 ^)' 7 0 = 7 ° 7 M , the Dirac equation for 
spinor T will be 

id^( 7 ^ ® I N ) + ^div( 7 ) <g> I N - 5 ^( 7 M ® /jv)v4“T a = . (5.4') 

Then multiplying (5.4) by 4 ~>T a from left and (5.4') by T 0 \T from right and adding the 
obtained equations, we get div(j a ) = div(j) = 0 if spinor T obeys the Dirac equation 
(5.4). 

The question now is how to connect the sought relativistic bound states with the sys¬ 
tem (5.4)—(5.5). To understand it let us apply to the experience related with QED. In the 
latter case lagrangian looks like (5.1) with changing group SU(-ZV)—>U(1) so T will be just 
a four-dimensional Dirac spinor. Then, as is known (see, e. g. Ref. 1231 ). when passing 
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on to the nonrelativistic limit the Dirac equation (5.4) converts into the Pauli equation and 
further, if neglecting the particle spin, into the Schrodinger equation, parameter /xq becom¬ 
ing the reduced mass of two-body system. The modulo square integrable solutions of the 
Schrodinger equation just describe bound states of a particle with mass /xo, or, that is equiv¬ 
alent, of the corresponding two-body system. Historically, however, everything was just 
vice versa. At first there appeared the Schrodinger equation, then the Pauli and Dirac ones 
and only then the QED lagrangian. In its turn, possibility of writing two-body Schrodinger 
equation on the whole owed to the fact that the corresponding two-body problem in classi¬ 
cal nonrelativistic (newtonian) mechanics was well posed and actually quantizing the latter 
gave two-body Schrodinger equation. Another matter was Dirac equation. Up to now no¬ 
body can say what two-body problem in classical relativistic (einsteinian) mechanics could 
correspond to Dirac equation. The fact is that the two-body problem in classical relativistic 
mechanics has so far no single-valued statement. Conventionally, therefore, Dirac equation 
in QED is treated as the relativistic wave equation describing one particle with spin one half 
in an external electromagnetic field. 

There is, however, one important exclusion - the hydrogen atom. When solving the 
Dirac equation here one considers mass parameter /xo to be equal to the electron mass and 
one gets the so-called Sommerfeld formula for hydrogen atom levels which passes on to 
the standard Schrodinger formula for hydrogen atom spectrum in nonrelativistic limit (for 
more details see, e. g. Ref. 1231 and also Subsection 6.2 ). But in the Schrodinger formula 
mass parameter /xo is equal to the reduced mass of electron and proton. As a consequence, 
it is tacitly supposed that in Dirac equation the mass parameter should be equal to the same 
reduced mass of electron and proton as in Schrodinger equation. Just the mentioned re¬ 
duced mass is approximately equal to that of electron but, exactly speaking, it is not the 
case. We remind that for the problem under discussion (hydrogen atom) the external field is 
the Coulomb electric one between electron and proton, essentially nonrelativistic object in 
the sense that it does not vanish in nonrelativistic limit at c —> oo. If now to place hydrogen 
atom in a magnetic field then obviously spectrum of bound states will also depend on pa¬ 
rameters decribing the magnetic field. The latter, however, is essentially relativistic object 
and vanishes at c —► oo because, as is well known, in the world with c = oo there exist no 
magnetic fields (see any elementary textbook on physics, e. g. Ref. EU). But it is clear that 
spectrum should as before depend of /xo as well and we can see that /xo is the same reduced 
mass as before since in nonrelativistic limit we again should come to the hydrogen atom 
spectrum with the reduced mass. So we can draw the conclusion that if an electromagnetic 
field is a combination of electric Coulomb field between two charged particles and some 
magnetic field (which may be generated by the particles themselves) then there are certain 
grounds to consider the given (quantum) two-body problem to be equivalent to the one of 
motion for one particle with usual reduced mass in the mentioned electromagnetic field. 
As a result, we can use the Dirac equation for finding possible relativistic bound states for 
such a particle implying that this is really some description of the corresponding two-body 
problem. 

Actually in QED the situaion is just as the described one but magnetic field is usually 
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weak and one may restrict oneself to some corrections from this field to the nonrelativistic 
Coulomb spectrum (e.g., in the Zeeman effect). If the magnetic field is strong then one 
should solve just Dirac equation in a nonperturbative way (see, e.g. Ref. (25|). The lat¬ 
ter situation seems to be natural in QCD where the corresponding magnetic (colour) field 
should be very strong because just it provides linear confinement of quarks as we shall see 
below (see also Section 9). 

At last, we should make an important point that in QED the mentioned electromagnetic 
field is by definition always a solution of the Maxwell equations so within QCD we should 
require the confining SU(3)-field to be a solution of Yang-Mills equations. Consequently, 
returning to the system (5.4)—(5.5), we can suggest to decribe relativistic bound states of 
two quarks (mesons) in QCD by the compatible solutions of the given system. To be more 
precise, the meson wave functions should be the nonperturbative modulo square integrable 
solutions of Dirac equation (5.4) (with the above reduced mass po) in the confining SU(3)- 
Yang-Mills field being a nonperturbative solution of (5.5). In general case, however, the 
analysis of (5.5) is difficult because of availability of the nonabelian current J of (5.6) in 
the right-hand side of (5.5) but we may use the circumstance that the corresponding modulo 
square integrable solutions of Dirac equation (5.4) might consist from the components of 
form T ; ~ r aj e with some ctj > 0, f3j > 0 which entails all the components of the 
current J to be modulo << 1 at each point of Minkowski space. The latter will allow us to 
put J « 0 and we shall come to the problem of finding the confining solutions for the Yang- 
Mills equations of (5.5) with J = 0 (according to Section 1 such solutions are supposed 
to be spherically symmetric and to contain only the components of the SU(A r )-field which 
are Coulomb-like or linear in r) and after inserting the found solutions into Dirac equation 
(5.4) we should require the corresponding solutions of Dirac equation to have the above 
necessary behaviour. Under the circumstances the problem becomes self-consistent and 
can be analyzable. 

It is clear that all the above considerations can be justified only by comparison with 
experimental data but now we obtained some intelligible programme of further activity. So 
let us pass on to its realization. 

5.2 Role of Diagonal Gauge 

Let us in detail write out .4 /; = of the Dirac operator from (5.2) employing the 

SU(iV)-Lie algebra realizations from Appendix B. We obtain at N = 2, 3,4 respectively 
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/ + A„ + A„ 


z 1 


_ /l3 1 46 _ 4 15 


V 


-22 

23 \ 

-24 

25 

46 415 

26 

* 

^6 

- Al-Al + A 15 ) 


(5.7) 


M 

2 ? 

^2 ^4 

^3 ^5 

with 24 = ^ + yl® — z(Aj; + d 42 ), ^2 = ^4 4 + A^f - i(A 3 + d 44 ), Z 3 = A 7 ^ — A 44 - 
+ A™), z 4 = ^ + 7L 44 — *(^ — A^ 0 ), 25 = A 4 — A ^ 3 — i(A 3 — A 44 ), Z 6 = 
A 4 — A® — i(A 2 — A 42 ), where (*) signifies complex conjugation. Then it is natural to 
put A“ = 0 with a = 1,2 at At = 2, A“ = 0 with a = 1, 2 ,4, 5, 6 , 7 at N = 3, A“ = 0 
with a = 1, 2,4, 5, 7, 8,9,10,11,12,13,14 at N = 4 so long as the Dirac equation (5.4) in 
such a gauge takes the simplest form. We further call this gauge diagonal one. Really Dirac 
equation (5.4) in diagonal gauge splits into the system of Dirac equations for components 
. Namely, at N = 2 the system is 


il e K 


Of, - - igAl 


il e K 


~ 2 ^b7 a 7 + ig^l 


*1 = 00*1 , 


*2 = 00*2 , 


(5.8) 


while at N = 3 it is 

*7 e E» 


d M ~ -uj fi abl a l b -iglA^A -4=4 


1 


9 >„ - \^!xabl a l b - ig (-Al + =Aj 
- ig 


and, at last, at N = 4 the corresponding system is 

il e E» 


*1 = 00*1 , 

*2 = 00*2 , 
*3 = , 


(5.9) 


d, - \^abi a i b - ig (4 3 + Al + A 15 ) 


*1 = 00*1 , 


il e E» [3 M - \u lxabl a l b - ig [~Al + A% - A 1 ?) 
H e E £ [a M - \^abl a l b ~ ig { Al - Al - Al 5 ) 
d, ~ - ig {~Al -Al + A 15 ) 


il e E» 


*2 = 00*2 , 
*3 = 00*3 , 
*4 = 00*4 • 


(5.10) 


It is cleai - that in diagonal gauge the SU (IV)-Yang-Mills fields arc described by matrices 
A /; = A“T a taking their values in the Cartan subalgebra of the conforming SU()V)-Lie 
algebra (see Appendix B). 

Due to that da* A da* = —dx u A da*, we have A A A = A“A* [T a , T h \ da* A da*, 

fi < v. Consequently for the S U (A r ) - Yang-M i 11 s with values in the Cartan subalgebra 
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A A A = 0 since commutator [T a , 7),] for diagonal matrices is always equal to zero and the 
Cartan subalgebras of SU(A r )-groups just consist from diagonal matrices (see Appendix 
B). Accordingly the curvature matrix (field strength) F = dA + gA A A = dA while the 
right-hand side of the Yang-Mills equations (5.5) (at J = 0) is identically equal to zero 
since matrix *F is also diagonal and then *F A A = A A *F. This fact strongly simplifies 
the task of searching for confining solutions because the equations (5.5) convert into 

d*F = 0. (5.11) 

It should be emphasized that such a simplification is dictated by the wish to obtain the 
simplest form for the Dirac equation (5.4). Clearly, all the results obtained in diagonal gauge 
will hold true in any gauge connected with the diagonal one by some gauge transfomation 
due to the fact of gauge invariance of the Yang-Mills equations (5.5). It turns out, however, 
that the confining solutions obtained in this way has the property of uniqueness in a certain 
sense which will be discussed in Section 8 and now let us pass on to finding confining 
solutions. 

5.3 U(l)-case 

We have already discussed this case in Section 2, where the corresponding confining solu¬ 
tion was found in the form (2.6). But it could seem that when searching for the solution 
the ansatz used was not the most general one. Really, we took the ansatz in the form 
A = A t {r)dt + A v (r)dp. It seems that the most general form is A = A t (r)dt + A v (r)dr + 
A.#(r)dd + A<p(r)d<p. Let us discuss it in more details. 

For the latter ansatz we have F = dA = —d r A t dt A dr + d r A^dr A dd + d r A ip dr A dip 
for an arbitrary A r (r). Then, according to (A.6), we obtain 

*F = ( r 2 sin &)d r Aid'd A dp + sin dd r A.^dt A dp - — r rd r Au,dt A (5.12) 

sm v 

which entails 

d* F = sin dd r (r 2 d r A t )dr A dd A dp — sinr 9E%A#dt A dr A dp- 

cos i!)d r Addt A did A dp H- —-d^A^dt A dr A did = 0 , (5.13) 

sm v 

wherefrom one can conclude that 

d r (r 2 d r A t ) = 0, d 2 A v = 0 , (5-14) 

d 2 r A d = drA# = 0 . (5.15) 

This yields the solutions (2.6) while we draw the conclusion that A& = C\ with some 
constant Cj. But then the Lorentz condition (2.2) for the given ansatz entails 

sinfh9 r (r 2 A.) + 9^(sinr?A^) = 0, 
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or 

d r (r 2 A r ) + cot dA$ = 0, (5.16) 

which entails A r = — A-q cot 0/r + C' 2 /r 2 with a constant C' 2 . But the confining solutions, 
as we accept in the given paper, should be spherically symmetric and contain only the 
components which arc Coulomb-like or linear in r, so one should put C\ = C '2 — 0. 
Consequently, the ansatz A = A t (r)dt + A v {r)dip is most general and we can consider 
A r = ,4,y = 0 without loss of generality. 

Let us describe one class of the confining nonspherically symmetric solutions of the 
Maxwell equations (2.3) (at J = 0) that can be obtained with the aid of the ansatz A = 
A t (r)dt + A^(r, d)dp, i. e. now we consider the component A, p depending also on d. 
It is evident that the Lorentz condition (2.2) is automatically fulfilled for the given ansatz. 
We shall have F = dA = —d r A t dt A dr + d r A v dr A dip + d^A^dd A dip and *F = 
r 2 sin dd r A t dd A dip - F^d r A v dt Add + ^F^d^A^dt A dr. Then Eq. (2.3) (at J = 0) 
entails 

d r (r 2 d r A t ) = 0, (5.17) 

+ sin dd^ (= 0 . (5.18) 

\ sin d ) 

We shall not here discuss the general form of the solution for Eq. (5.18) and only write out 
the possible solution of (5.17)—(5.18) which is useful to us in the present paper in the form 
slightly modifying (2.6) 


A t = —b A , A,„ = hr + B — K cos d (5.19) 

r 

with some constants a, b, A, B, K parametrizing solution. 

5.4 N = 2 

Remarks done in previous subsection about the Lorentz condition will hold true for any 
group SU(IV) (see Subsection 8.1), so at N = 2 we put A'F } = 0. After this we search for 
the solution of (B. 3) (at J = 0) in the form A = A t (r)dt + A v {r)dip with A t = Af 
Along the above lines it is then easy to come to the system 

d r {r 2 dr At) = 0, d 2 r A v = 0 , (5.20) 

and we write down the solutions of (5.20) needed to insert into (5.8) 

A 3 t =-- + A, (5.21) 

r 

A\^ = br + B 

with some constants a, A, b, B parametrizing solutions. 


(5.22) 
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Class of the confining nonspherically symmetric solutions of (7i.3j can be obtained with 
the aid of the ansatz A = A t {r)dt + A v (r. d)d<£>. We shall have 

d r {r 2 d r A t ) = 0, (5.23) 

r 2 d 2 Ap + sin dd# (d^A^) = 0 . (5.24) 

\ smv J 

It is clear that Eq. (5.23) gives the same solution of (5.21) while one possible solution of 
(5.24) useful to us in the present paper is 

A 3 = — K cos d + br + B (5.25) 

with some real constants K. b , B parametrizing solution. 


5.5 N = 3 


In the given case we put A 8 ' 8 = 0. After this the ansatz in the form A = A t [r)dt+Ap(r)dtp 
with A t jp = .4A. 3 +.4* , A 8 yields (at J = 0) the solutions (2.13). The corresponding class 
of the confining nonspherically symmetric solutions of (B. 3) (at J = 0) can be obtained 
with the aid of the ansatz A = A t (r)dt + A v (r, id)d(p which gives the same component A t 
as in (2.13) while 


4 3 -I- — A 8 - 

* Vs v ' 

- 4 3 4 - —A 8 

* v/3* 


K\ cosd + b\r + B\ , 
-K 2 cos d 4 -b 2 r + B 2 , 



{K\ 4- K 2 ) cos d - ( 6 i 4- b 2 )r 


[Bx + B 2 ) 


(5.26) 


with some real constants Kj , bj , Bj parametrizing solution. 


5.6 N = 4 

We put A 8 ^' 15 = 0 and at J = 0 the ansatz A = A t (r)dt + A v (r)d(p with A t)ip = 
Al .74 4 - AfpTg 4 - .4) (47) 5 gives rise to the solutions of ( BA) in the form 


AS + Af 


+ A ] 5 = —— 4- Ai ,-Af + A 8 - AY = —— 4- A 


a i 


15 


02 


At — A® — Aj 5 = —+ A 3 , —A 8 — Af + A } 5 = — + + ° 3 — {A\ 4- A 2 + A 3 ) , 

A% + A® 4 - A = b\r 4- B\ , — A 3 4 - A — A ^ 5 = b 2 r A B 2 , 

Ap — A®—Ap 5 = 63 r4-^3) — Ap — A® 4-A * 5 = — {bi+b 2 +b 3 )r—(Bi+B 2 +B 3 ), (5.27) 

where real constants ctj , Aj, bj, Bj parametrize the solutions. The corresponding class of 
the confining nonspherically symmetric solutions of (B. 3) (at J = 0) can be obtained with 
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the aid of the ansatz A = A t (r)dt + A v (r, d)dp which gives the same component A t as in 
(5.27) while 

A 3 v + A^ + A 25 = b\r + Bi - Ki cos d , -A^ + A® - A * 5 = b 2 r + B 2 - iT 2 cos d , 
A % - A % - A !p = hr + B 3 - K 3 cos d , 

-4 - K + 4 5 = -(h + &2 + & 3 )r - (Bi + B 2 + Bi) + (K\ + K 2 + A 3 ) cos d (5.28) 
with real constants Kj. bj, Bj. 


6 Spectrum of Bound States in the Coulomb-Like Case 

The question now is how to find the modulo square integrable solutions of Dirac equa¬ 
tion (5.4) when inserting the confining solutions described in previous Section into it. We 
shall need some results about spectrum of the euclidean Dirac operator on the unit two- 
dimensional sphere § 2 in the form obtained in Refs. ESI- 

6.1 Results from the Black Hole Theory about Eigenspinors of the (Twisted) 
Euclidean Dirac Operator on S 2 

When separating variables in (5.4) (see next Subsection) there naturally arises the euclidean 
Dirac operator Vq on the unit two-dimensional sphere § 2 and we should know its eigenval¬ 
ues with the corresponding eigenspinors. Such a problem also arises in the black hole theory 
while describing the so-called twisted spinors on Schwarzschild and Reissner-Nordstrom 
black holes and it was analysed in Refs. EE so we can use the results obtained therein 
for our aims. Let us adduce the necessary relations. 

Let us consider 2A:-dimcnsional (pseudo)riemannian manifold M for which // 1 (A/, Z 2 ), 
the first cohomology group with coefficients in Z 2 , is equal to zero while H 2 [M. Z), the 
second cohomology group with coefficients in Z, is equal to Z. Then standard topologi¬ 
cal results 118j J9 [20j| say to us that over M there exists the only so-called Spin-structure 
whereas there is countable number of complex line bundles over M. As a consequence, 
each complex line bundle can be characterized by an integer n which in what follows will 
be called its Chern number. Under this situation, if denoting S(M ) the only standard spinor 
bundle over M and £ the complex line bundle with Chern number n, we can construct ten- 
sorial product S(M ) <g) £. Under the circumstances we obtain the twisted Dirac operator 
V n : S(M) 0 £ S(M ) <8> £, so the eigenvalue equation for corresponding spinors 4> as 
sections of the bundle S(M ) 0 £ may look as follows 

T>n& = A$, (6.1) 

and we can call (standard) spinors corresponding to n = 0 (trivial complex line bundle £) 
untwisted while the rest of the spinors with n/0 should be referred to as twisted. 



Intersection of Black Hole Theory and Quantum Chromodynamics... 


23 


From general considerations mmm the explicit form of the operator V n in local 
coordinates x 1 ' on a 2/.'-dimensional (pseudo)riemannian manifold can be written as follows 

V n = = iYE^dfj, - ^o; M a 67 a 7 b - ieA^), a < b, (6.2) 

where A = A^dx 11 is a connection in the bundle £ and the forms uj a f, = oj^dx^ obey the 
Cartan structure equations de a = u a b A e b with exterior derivative d, while the orthonormal 
basis e a = e^dxJ' in cotangent bundle and dual basis E a = E£d^ in tangent bundle are 
connected by the relations e a (E} ) ) = At last, matrices 7 “ represent the Clifford algebra 
of the corresponding quadratic form in C 2k . Below we shall deal only with 2D euclidean 
case of the unit sphere S 2 (k = 1, quadratic form Q 2 = Xq + xf). 

As for the connection A jt in bundle £ then the suitable one can be found, for example, 
in Refs. (SJ and is 

Tl 

A = Andx ^ = -cos ddtp . (6.3) 

e 

Under the circumstances, as was shown in Refs. 0, integrating F = dA over the unit 
sphere S J gives rise to the Dirac charge quantization condition 

[ F = 47T— = 47T q (6.4) 

is 2 e 


with magnetic charge q, so we can identify the coupling constant e with electric charge. 

As was discussed in Refs. <261 . the natural form of V n in local coordinates '0. ip on the 
unit sphere S 2 looks as follows 


V n = —ia\ 


icr2 <9tf + io-3^—7 ( d v — \01cr 3 cos d + in cos d 
sm-d V 2 


+ ——rcri(73 d v -——<ti(72 + ino\U 3 cot d (6.5) 

sm d 2 

with the Pauli matrix aj (see Appendix B), so that rr 1 V n = —V n a\. As is not complicated 
to see, the operator V n has the form (6.2) with 7 0 = —%<7\ a 2 , 7 1 = —z<ti 173 , e° = dd, 
e 1 = sin ! dd(p,E 0 = 8$, E\ = d^/ sind, cjqi = cos dd<p, A^dx^ = —^ cos dd<p. 

The equation (6.1) was explored in Refs. EH. Spectrum of D n consists of the numbers 
A = ±yj(l + l) 2 — n 2 with multiplicity 2(1 + 1) of each one, where l = 0,1, 2,..., I > \n\. 
Let us introduce the number m such that — l < m < l + 1 and the corresponding number 
m' = m — 1/2 so \m'\ < l + 1/2. Then the conforming eigenspinors of operator V n are 


4 > = 


$1 

q > 2 


= = 


C 


' pk \ pk 

r m'n—1/2 ra'n+1/2 

pk -J- pk 

, r m!n—1/2 ' Jr m'n+1/2 , 


D —im cp 


( 6 . 6 ) 


with the coefficient C = 


These spinors form an orthonormal basis in L^S 2 ) for 
each n and are subject to the normalization condition 


7 T 27 T 


4F4>dD = J /(|^| 2 + |4> 2 | 2 )sin^ = l, 


0 0 


(6.7) 
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where (f) stands for hermitian conjugation. As to functions P^, n , (cost?) then they can be 
chosen by miscellaneous ways, for instance, as follows (see, e. g., Ref. EH) 


P m , n /(cosi?) = i 


E 


m -j -n 


{k — m')\(k — n')\ f 1 + cos■&\ 2 


(k + m')\(k + n')\ \1 — cost? 

(k + /l—cosilV 


j=max(m',n') ~ j)! (J ~ ™')! (j - n ')! 


( 6 . 8 ) 


with the orthogonality relation at m', n’ fixed 

7 T 

J P m'n‘ ( cos A P m'n' (cOS 1?) sin MV = EE_<W , 

0 

where (*) signifies complex conjugation. It should be noted that square of V n is 


(6.9) 


2 in cos dcL — n 2 




sin 2 i? 


1 2 

+ m • 2 - n 

sin 1 ? 


with 


2 cost? 1 

Do = —A s 2 + o~2 g~3 . Q ip + 


1 

i — r 2 n n 1 1 2 n 7 > 

sm 1 ? 4 sm d 4 


while laplacian on the unit sphere is 


A s 2 = sin $9$ + , 2 d 2 = dj + cot ^ + . \ > 

smt? snrt? ^ sm z i9 ^ 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 


so the relation (6.10) is a particular case of the so-called Weitzenbock-Lichnerowicz for¬ 
mulas (see Refs. [' 18. 19, 20]). Then from (6.1) it follows D 2< b = A 2 <b and, when using the 


ansatz 4> = P(i?)e 


—im cp _ 


Pi 

p? 




, P \-2 = Pi, 2 ($)> the equation D 2 <3? = A 2< b turns 


into 


— da — COt + 


m' 2 + n 2 + 4 — 2m!n cos il 777 / cos t? — 


sin 2 t? 


A 2 + n 2 - 


+ 


n 


sin 2 i) 


CJl\P = 


P 


(6.13) 


wherefrom all the above results concerning spectrum of V n can be derived 626]. 

6.2 U(l)-case 

We should insert the confining solutions (2.6) into Dirac equation (5.4), where T will be 
just a four-dimensional Dirac spinor, and let us employ the ansatz 


4 / = e lujt r 1 


Pi(r)4>(i?,<p) 

F 2 (r)ai^('&,(p) 


(6.14) 



















Intersection of Black Hole Theory and Quantum Chromodynamics... 


25 


/<f>, \ 

with a 2D spinor <f> = I I. Then, after a simple matrix algebra computation, we can get 
from (5.4) the system 



+ -V 0 - 
r 


°2 

sin v 



*(A*o 


c) 


-F 2 $, 



+ -22o 
r 


°2 

smv 



—F2& 1$ = 
r 


-i(go + c)-Fiai& 

r 


(6.15) 


with c = lo — g(—a/r + A) while the euclidean Dirac operator Vq on the unit sphere § 2 is 
given by (6.5) at n = 0. It is not complicated to check that at 6 / 0, B / 0 the variables 
r and v) are not separated. Under this situation we shall at first restrict ourselves to the case 
6 = B = 0 since under the circumstances we can solve Eq. (5.4) exactly. Really we employ 
the ansatz (6.14) and obtain the system (due to the fact that a{D 0 = —Vq(t{) 



i(Ho - c)F 2 , 


(d r - F 2 = -i(m ) + c)Fi (6.16) 

with an eigenvalue A for the eigenspinor <I> of the above operator Vq, X = ±(l + 1) e 
Z\{0} , l = 0,1, 2... (see previous Subsection). 

Let us now employ the ansatz 

Fi = yjn 0 - (w - gA ) r a e~ 0r [fi{x) + / 2 (x)], 

F 2 = i\JA 0 + (w - gA) r a e~ 0r [fi{x) - f 2 (x)] (6.17) 

with a = \/X 2 — g 2 a 2 , (3 = \J — (cu — gA) 2 , x = 2/3r. 

Then, inserting the ansatz into (6.16), adding and subtracting equations give rise to 


flxf[ + Yfi + Zf 2 = 0 , (6.18a) 

f3xf 2 - (ixf 2 + lo/ 2 + Z {) }\ = 0 , (6.186) 

where prime signifies the differentiation with respect to x, Y, Yq = a/3 ga{uj — gA), 
Z , Zo = A/3 ± gag o- From (6.18), if using the relations YY$ — ZZq = 0, Y + Yo = 2a/3, 
one yields the second order equations in x 

x fi + (! + 2a - x)f[ - ^/i = 0 , (6.19a) 

x} 2 + (1 + 2a - x)f 2 - (l + Fj f 2 = 0 , (6.196) 

that are the Kummer equations (confluent hypergeometric equations in another terminology, 
see, e. g. Ref. t29l ) and for (6.19a) the only finite solution at 0 and at infinity not strongly 
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increasing is the Laguerre polynomial L£(x) with n 
spectrum 


ui = gA ± g 0 


1 + 


S ' 2 a 2 


(n + \/A 2 - 


-17/3 = 0,1,2, 


-1/2 


g 2 a 2 ) 2 


This gives the 


( 6 . 20 ) 


wherefrom it is clear that constant A only shift the origin of count for energy and we can 
consider A = 0. Further, putting /i = CL 2 ' y (x) with some constant C, from (6.18a) at 
n > 0 we find 


h 



2a+l 
n —1 


(x)-YLl a (x) 


because [L^(x)]' = —L 2 7V (x) ll28l . that entails 


Y 


F\ = Cy/g 0 -ur a e 0r (1 — — ) L 2 n a (x ) + FxL 2 v a F (x) 


7 


z 


F 2 = iC^+ZJr a e~ 0r 


Y 


Z 


(3 


1 + 17 L l a ( x ) - ^Ll a _\\x) 


( 6 . 21 ) 


The case n = 0 should be considered separately. We here have Y = 0, Yq = 2a/3 = 2gau, 
u = ±go\/\ 2 — g 2 a 2 / |A|, f\ = CLq°(x) = C. Further ZZo = 0 = (A/3) 2 — (gag o) 2 
which entails |A|/3 = g\a\go. Then at a > 0, A > 0 we get Z = A/3 + gag o = |A|/3 + 
g\a\g 0 = 2gag 0 > 0, Z 0 = A/3 - gag 0 = |A|/3 - g\a\g 0 = 0, f 2 = -CY/Z = 
—CZq/Yq = 0 and 


Fi = Cy/go - w r a e~ 0r , F 2 = iCy/g 0 + u r a e~ 0r . (6.22) 


At a > 0, A < 0 we obtain Z = X/3+gago = —\\\(3+g\a\go = 0 > 0, Zo = —2gago < 0, 
h = —CZq/Yq = Cg 0 /u = C \\\/(±^/\ 2 - g 2 a 2 ) and 

F, = CVg^r a e~ 0r (l + ^) ,F 2 = iC^+Z3r a e~ 0r (l - . (6.23) 

At a < 0, A < 0 we get Z = A/3 + gago = — |A|/3 — g\a\go = — 2g\a\go < 0, Zo = A/3 — 
gago = — |A|/3 + g\a\go = 0, f 2 = 0 and F ])2 are given by (6.22). At last, at a < 0, A > 0 
we have Z = A/3+firap 0 = |A|/3-p|a|/r 0 = 0, Z 0 = \{3-gag 0 = |A|/3+fir|a|p 0 = 2g\a\go, 
f 2 = —CZq/Yq = Cgo/u and Fi ) 2 are given by (6.23). 

To describe relativistic bound states we should require dr e L 2 (M 3 ) at any t £ R and 
one can accept the normalization condition for F \, F 2 in the form 

roc 

/ (|Fi| 2 + |F 2 | 2 )dr = l (6.24) 

Jo 

with taking into account the condition (6.7) so that C can be determined from the relation 
(6.24). 

It is the expression (6.20) that is in essence the Sommerfeld formula mentioned in Sub¬ 
section 5.1. It should be noted, however, that standard parametrization in the Sommerfeld 
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formula adduced in all the monographs (see, e. g., C3E3I2B) is somewhat different 
from that in formula (6.20). The fact is that the standard approach uses the orthonor¬ 
mal basis in L 2 (§ 2 ) different from the basis of the eigenspinors of the Dirac operator Vq. 
Namely, one uses the eigenbasis of the operator K = aL + 1, where a = (o \. a 2 , a 3 ) and 
L = — i(f X d/dr) is the angular momentum operator. The eigenvalues of operator K are 
numbers of the form A 7 = ±(j + 1 / 2 ) with j = 1 / 2 ,3/2,..., so j defines the total angular 
momentum J = L + S with spin momentum S = a/2, that is J 2 has the eigenvalues 
j ( j + 1). In standard approach just A' stands in (6.20) instead of A. It is evident that both 
formulas (with A or A 7 ) reproduce the same spectrum but the corresponding wave functions 
will be slightly different in angular part depending on $, ip since the Dirac operator 'D () does 
not commute with operator K. There are at least two reasons why historically operator K 
was employed rather than Vq. The first one is that in standard approach one solves Dirac 
equation (5.4) in Cartesian coordinates without transition to spherical ones and under this 
situation there naturally arises just operator K whereas when passing on to the spherical 
coordinates there naturally would arise just operator Vq while separating variables. The 
second reason is that the whole formula (6.20) is not necessary in nonrelativistic systems 
such as hydrogen atom, positronium and so on. It is enough to restrict themselves to a 
few terms of expansion in g 2 for (6.20) to obtain collections (fine structure) to the non rela¬ 
tivistic spectrum which arc in concordance with experiment and one may use the notion of 
spin as essentially nonrelativistic phenomenon. Another matter are quarks in mesons that 
should probably be considered essentially relativistic objects and it is the most natural to 
write down Dirac equation (5.4) in Minkowski spacetime with spherical coordinates in its 
spatial part. It is the latter way that we pursue in present paper. 

At the end of this Subsection we can slightly generalize the results obtained if inserting 
solution (5.19) at b = B = 0 into Dirac equation (5.4). Then it is not complicated to see we 
shall get the similar relations providing that K = k/g with k E Z, i. e., k is integer number. 
But now we should consider the spinor <I> of (6.14) to be the eigenspinor of the twisted 
euclidean Dirac operator on the unit sphere § 2 with the Chern number k (see Subsection 
6.1) and the eigenvalues A should be, accordingly, replaced by A = ± \J(l + l ) 2 — A: 2 , 
l > |/e|. Physically the corresponding configurations describe the Dirac monopole ones 
with magnetic charge P = k/g so the corresponding wave functions should be modified in 
obvious way. 


6.3 N = 2 


After inserting the solution (5.21) (with A = 0) into (5.8) for both the equations we employ 
the ansatz 


T, = 


— 1 


i4(r)<4(44 \ 

F j2 {r)ai^j(^,(p) ) ’ 


with a 2D spinor 4 


4i 

4-2 


which entails the systems 


3 = 1,2 


(6.25) 
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^<9 r - F j2 = ~i(no + Cj)Fj i (6.26) 

with an eigenvalue A j for the eigenspinor <l>j of the Dirac operator T>q, \j = ±(lj + 1) £ 
Z\{0} , lj = 0,1,2... Besides 

ci=ui + ga/r,c 2 =U 2 - ga/r , (6.27) 


so that the energy spectrum uj of particle is given by the relation lo = ui\ + ui 2 . Acting along 
the same lines as in previous Subsection we obtain the spectrum of particle in the form 




-1/2 


- = ± 

1 1 9 a 

± 

1 9 a 

Po 

(ni +J\\-g 2 a 2 ) 2 _ 


(n 2 + JX 2 - g 2 a 2 ) 2 _ 


-1/2 


(6.28) 


where the number m ;2 = 0 , 1 , 2 ,.... 

If K / 0 in (5.25) at b = B = 0 then when inserting (5.21), (5.25) into (5.8) we shall 
get the si mi lar spectrum (6.28) providing that K = k/g with k £ Z, i. e., k is an integer 
number. But now we should consider the spinor <f>j of (6.25) to be the eigenspinor 4/ ? of the 
twisted euclidean Dirac operators V±k on the unit sphere § 2 (see Subsection 6.1), respec¬ 
tively, with the Chern numbers ±k and the eigenvalues A j should be, accordingly, replaced 
by Ai = ±\/{li + l ) 2 — F 2 , h ^ |fc|> A 2 = ±^((2 + l ) 2 — k 2 , l 2 > \k\. Physically the 
corresponding configurations of SU(2)-field describe the Dirac-like monopole ones with 
magnetic charges, conformably, P\ = k/g, P 2 = —k/g but the total (nonabelian) magnetic 
charge of the given configurations remains equal to Pi + P 2 = 0. 

The corresponding wave functions are not complicated to be written out on the analogy 
of U(l)-case of previous Subsection but we shall not dwell upon it. It should be only noted 
that the condtion (6.24) should be replaced by 



1 

2 ’ 


3 = 1,2. 


6.4 N = 3 


We here use the solution for A t of (2.13) (with A\= A 2 = 0) and for all three equations in 
(5.9) we employ the ansatz 


'Pj = e iuJjt r 1 


Fjiir^j^y) \ 

Fj 2 (r)<n<bj(ti,<p)J 


,3 = 1 , 2,3 


with a 2D spinor ( I) ? 


( enta ^ s the systems 



Fj 1 = i(g 0 


c j)Fj2 1 


(6.29) 


d r — 




r 


F j2 = -i(no + Cj)Fji 


(6.30) 
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with an eigenvalue A j for the eigenspinor <J>j of the Dirac operator Vq, A j = ±(lj + 1) G 
Z\{0} , lj = 0,1,2... Besides 

ci = uji + gai/r, c 2 = to 2 + ga 2 /r, c 3 = w 3 - g(ai + a 2 )/r , (6.31) 


so that the energy spectrum c o of particle is given by the relation ui = uj\ + u > 2 + cc 3 which 
we obtain in the form 



r n 

-1/2 

- 

- = ± 

1 1 

± 

1 1 g2a2 

go 

(m + J Af - g 2 af) 2 _ 


(n 2 + J\l - g 2 al) 2 _ 


± 


1 + 


_ g 2 (ai + Q2) 2 _ 

( n 3 + \A| -g 2 (a 1 + o , 2 ) 2 ) 2 


-1/2 


(6.32) 


where the number ni ;2)3 = 0,1, 2,.... 

If Kj / 0 in (5.26) at bj = Bj = 0 then when inserting (2.13), (5.26) into (5.9) we shall 
get the similar spectrum (6.32) providing that Kj = kj/g with kj G Z, i. e., kj arc integers. 
But now we should consider the spinor 4>j of (6.29) to be the eigenspinor <h ; of the twisted 
euclidean Dirac operator on the unit sphere S 2 (see Subsection 6.1), respectively, with 
the Chern numbers k = k 1. /r 2 , —{k\ + k 2 ) and the eigenvalues A j should be, accordingly, 

replaced by Ai = ±yj(l 1 + l) 2 - kf, h > \h\, A 2 = ±^J(l 2 + l) 2 - k%, h > \k 2 \, 
A 3 = i \J (/ 3 + l) 2 — (fei + fe 2 ) 2 , ( 3 > |fci + /c 2 |. The coiTesponding configurations of 
gluonic field describe the Dirac-like monopole ones with magnetic charges, conformably, 
P\ = ki/g, P 2 = k 2 /g, P 3 = — (k\ + k 2 )/g, but the total (nonabelian) magnetic charge of 
the given configurations remains equal to P\ + P 2 + P 3 = 0. 

The corresponding wave functions arc again not complicated to be written out on the 
analogy of U(l)-case of Subsection 6.2 but we do not dwell upon it and the condtion (6.24) 
should be replaced by 


roo 1 

/ (\F jl \ 2 + \F j2 \ 2 )dr = -,j = l,2,3. 

Jo 3 


6.5 N = 4 


If inserting the solution for A t of (5.27) (with ,4i = .4 2 = ,4 3 = 0) into (5.10) then for all 
four equations in (5.10) the ansatz of form (6.29) with j = 1,2, 3,4 will lead to the systems 
of form (6.30) with an eigenvalue A j for the eigenspinor <I> y of the Dirac operator Vq, where 
A j = ±(lj + 1) G Z\{0} , lj = 0,1, 2..., while 

ci = ix>i+gai/r,c 2 = u 2 +ga 2 /r,c 3 = u} 3 +ga 3 /r,c 4 = w 4 -g(ai+a 2 +a 3 ))/r, (6.33) 


so that the energy spectrum ui of particle is given by the relation u = + u> 2 + w 3 + 0 J 4 

and is obtained in the form 


u> 

go 


1 + 


g 2 a\ 


(m + yj A? - g 2 a\f 


-1/2 


1 + 


_ g^a |_ 

(n 2 + y/\l ~ g 2 a ? 2 ) 2 


1/2 
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-| -1/2 r -i -1/2 

± , 9 2 a 2 3 ± 1 | g 2 (ai + a 2 + a 3 ) 2 

(«3 + - fi r2a 'i) 2 _ _ (n 4 + -5 2 («i +«2 + a 3 ) 2 ) 2 _ 

(6.34) 

where the number 71^2,3,4 = 0,1, 2,.... 

If Kj / 0 in (5.28) at bj = Bj = 0 then when inserting (5.27), (5.28) into (5.10) 
we shall get the similar spectrum (6.34) providing that Kj = kj/g with kj G Z, i. e., 
kj arc integers, so the spinor 4>j is the eigenspinor ( I> ; of the twisted euclidean Dirac op¬ 
erator Vj, on the unit sphere § 2 (see Subsection 6.1), respectively, with the Chern num¬ 
bers k = k\. k‘ 2 , A’s, — (k 1 + h '2 + hi) and the eigenvalues A j should be, accordingly, re¬ 
placed by Ai = ±y/(h + l) 2 - kj, h > |Aq|, A 2 = ±\J(h + l) 2 - h > \k 2 \, A 3 = 

+ l) 2 — &3> ^3 > |^31 ^4 = 4= v/p4 + l) 2 — {k\ + k 2 + A/?) 2 , I 4 > |Aq + fc 2 + 

The coiTesponding configurations of SU(4)-field describe the Dirac-like monopole ones 
with magnetic charges, conformably, P\ = k\/g, P 2 = k 2 /g, P3 = k%/g, P 4 = — {k\ + 
A' 2 + hi)/g but as before the total (nonabelian) magnetic charge of the given configurations 
remains equal to P\ + Ik + 77 + A 4 = 0. 

The corresponding wave functions arc again not complicated to be written out on the 
analogy of U(l)-case of Subsection 6.2 but we do not dwell upon it while the condtion 
(6.24) should be replaced by 

roo \ 

J o (\ F ji\ 2 + \ F j 2 \ 2 )dr = -,j = 1,2,3,4. 

7 Spectrum of Bound States in the Coulomb-Linear Case 

7.1 U(l)-case 

We now should return to the system (6.15) at b 7A 0, B yA 0, A = 0. One may hope to 
obtain the almost exact solution of (6.15) if considering rr 2 <b ~ .sin 1 ) <I>. As follows from 
the estimate in Appendix D this condition is rather good fulfiled and when doing so we 
probably make an error retaining eigenvalues A of the euclidean Dirac operator Vq on the 
unit sphere S 2 instead of the eigenvalues of a less symmetric operator on § 2 whose form is 
unknown explicitly. 

Having accepted the mentioned condition we come to the system (owing to the fact that 
ciA) = —72 0 <ti, ct 2 cti = -cricr 2 ) 

A / B \ 

9 r + - ~ 3 (b + — j Fi= i(p 0 - c)F 2 , 

A / B\ 

dr - \-g(b-\ -) F 2 = -iifio + c)Fi (7.1) 

r V r / 


with c = (j + ga/r. 
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Now we employ the ansatz 

Fi = Pr a e~ p, '[fi{x) + f 2 (x)\,F 2 = iQr a e~ 0 r [fi(x) - f 2 (x)] (7.1') 

with a = \/{X — gB) 2 — g 2 a 2 , (3 = \J — w 2 + g 2 b 2 , P = gb + /?, Q = /io + oj, 
re = 2f3r. 

After this, inserting the ansatz into (7.1), adding and subtracting equations entail 


x-PQ/l + y/i + Zh = 0 , (7.2a) 

xPQf 2 - xPQh + Yof -2 + (z 0 - jPQx) fi = 0 , ( 7.2 b) 

where prime signifies the differentiation with respect to x, Y, Yq = [a/3 T gaiv + gab\Q ± 
g 2 abP , Z, Zq = [(A — gB)P ± papo)]/? 3= g 2 abP and YYq — ZZ§ = 0. 

From (7.2a)-(7.26) one yields the second order equations in x 


+ (1 + 2a - x)f[ + nfi = 0 , 


, Z 0 - b 0 x 
with &o = gbPQ/(3, k = PQ/Y and 


+ 1/2 = 0 


n = 


gbZ - (3Y 


PPQ 


which entails the equation for spectrum of i o 


(7.3) 

(7.4) 


(7.5) 


[ g 2 a 2 + (n + a) 2 ]uj 2 + 2(A — gB)g 2 abu>+ 

[(A - gB ) 2 ~(n + a) 2 ]g 2 b 2 - ^(n + a ) 2 = 0 , (7.6) 

that yields 

_ -(A -gB)g 2 ab±\[X 

g 2 a 2 + (n + a) 2 1 ’ 

with X = (X — gB) 2 g 4 a 2 b 2 — [g 2 a 2 + (n + a) 2 ]{[(X —gB ) 2 — (n + a) 2 ]g 2 b 2 — ^(n + a) 2 } 
and it is clear that the expression (7.7) passes on to (6.20) at b = B = 0 while (7.7) can be 
rewritten in a more symmetrical form 


uj = aj(n, l , A) 


—A g 2 ab ± (n + a)^J (n 2 + 2na + A 2 )^q + g 2 b 2 (n 2 + 2na) 
n 2 + 2roa + A 2 


(7.8) 


with A = A — gB = +(/ + 1) — pH. 

It is clear that according to (7.3) (which is the confluent hypergeometric equation) we 
should choose f\ = CL 2a (x) with the Laguerre polynomial L 2 q (x) and some constant C 
if n = 0,1, 2.... Then at n > 0 from (7.2a) we find 


h 


C r 
~Z - 


xPQL 2 n a + 1 (x)-YL 2 n a (x ) 


2a I 
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because [L^ix)]' = ESI , that entails 


Fl = CPr a e-f* [(l - L 2 n a (x) + ^-xL 2 n a _+ *(x) 


F 2 = iCQr a e 


a -pr 


1 + 


Y 

~Z 


Ll a (x) - ^xL't\Ux) 


At n = 0 we have f\ = CLq(x) = C, f 2 = —CY/Z, wherefrom 


(7.9) 


Fi = CPr a e~P r (l - ^ = CPr a e~ 0r (^ ~ j) , 

F 2 = iCQ r a e~P r (l + |) = iCQ (l + ^) , (7-10) 

inasmuch as = /JY at n = 0 according to (7.5). 

Hence we can see that T of (6.14) e L^R 3 ) at any t £ R and, conformably, 4' may 
describe relativistic bound states with the energy spectrum (7.8). 

Also it should be noted that the influence of the Dirac monopole configurations for 
U(l)-field when K ^ 0 in (5.19) can be treated by the same manner as in Subsection 6.2 if 
taking cr 2 <l> « sin /M> for the eigenspinor $ of the twisted euclidean Dirac operator T>i~ on 
the unit sphere S 2 with the conforming Chern number k. 

At last, constant C from (7.9)—(7.10) is defined by the condition (6.24) as before. 


7.2 N = 2 

The corresponding modifications of U(l)-case here are obvious so we shall only briefly 
describe them. The spectrum is given by ui = oj\ + w 2 with 

= 0Ji(n\, li, Ai) = 


—Ai g 2 ab ± (n\ + ai)\J{n\ + 2n 1 o: 1 + A\)$ + g 2 b 2 (n\ + 2n 1 a 1 ) 


n- 


T 2ni«i -T A\ 


x>2 — w 2 (n 2 , l 2l A 2 ) — 

-A 2 g 2 ab ± (n 2 + « 2 )^/(n 2 + 2 n 2 a 2 + A|)/r^ + g 2 b 2 (n 2 + 2 n 2 a 2 ) 


n 


9 + 2 n 2 a 2 + Ao 


(7.11) 


(7.12) 


with Ai = X\-gB = ±(li + l)-gB, A 2 = X 2 +gB = ±(l 2 +l)+gB, a\ = JA. 2 - g 2 a 2 , 


a 2 = \JA 2 — g 2 a 2 . The corresponding radial parts of wave functions (6.25) are given at 
rij = 0 by 

Fj 1 = CjP^e-^ (l ~y)’ F * = (l + , 


(7.13) 
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while at rij > 0 by 


Fj\ — CjPjr 


a Pi 



PjQj 

(Tj) H— --r^-r 


to) 




Fj2 = i( \jQ jr" j < i ‘ r 


1 + Y ) (Q) - ^ 2 - r j L l a j J -i\ r j) 

\ 3 ) 3 


(7.14) 


with the Laguerre polynomials L^.ijrj), rj = 2 for, where /3j = yMo — cu 2 + g 2 b 2 , P\ = 
gb + (3i, P 2 = -gb + fa, Qj = go + Y\ = («i/3i ~ gau t + gai&)Qi + g 2 abPi, 
Y 2 = (a 2 P 2 + gau 2 - ga 2 b)Q 2 + g 2 abP 2 , Z x = [(Ai - gB)Pi + gag 0 ]Qi + g 2 abPi, 
Z 2 = [{X 2 +gB)P 2 —gago\Q 2 +g 2 abp 2 - Also it should be noted that the quantum numbers 
rij are defined by the relations 


_ gbZ 1 - fi\ Y\ _ gbZ 2 + (5 2 Y 2 
Ul P 1 P 1 Q 1 ,n2 P 2 P 2 Q 2 

Further, Cj is determined from the normalization condition 

roo 1 

(l^jr'tl 2 + \Fj 2 \ 2 )dr = - . 


(7.15) 


(7.16) 


Consequently, we shall gain that in (6.25) T'j £ L))® 3 ) at any t. 6 R. and, as a result, 
T = (Ti. T 2 ) may describe relativistic bound states with the energy spectrum (7.11)- 
(7.12). 

Finally, it should be noted that the influence of the Dirac-like monopole configurations 
for SU(2)-field when K h 0 in (5.25) can be treated by the same manner as in Subsection 
6.3 if taking a 2 ^j ~ sin for the eigenspinor 4> ; of the twisted euclidean Dirac operator 
T>±k on the unit sphere § 2 with the conforming Chern numbers 


7.3 N = 3 


Let us adduce the coresponding results without going into details that can be easily recon¬ 
structed along the lines of Subsection 7.1. The spectrum is given by uo = u>i + u> 2 + W3 
with 

mi = = 

-Aig 2 aibi ± ( n\ + a\)^J ( n\ + 2n\ai + A\)g% + g 2 b\(n{ + 2n\ai) 

n 2 + 2n\a\ + A 2 

LO 2 = <-02(n>2, h-, A 2 ) = 

-A 2 # 2 a 2 6 2 ± (n 2 + a 2 )\J ( n\ + 2n 2 a 2 + A|)/x| + g 2 b\{n\ + 2n 2 a 2 ) 

r?| + 2n 2 a 2 + 

^3 = <^ 3 ( 77 . 3 , ( 3 , A 3 ) = 


(7.18) 
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-A 3 g 2 a 3 b 3 ± (n 3 + a 3 )^(raf + 2 n 3 a 3 + A§j/xg + 5 2 6 §(n§ + 2ra 3 a 3 ) 

+ 2n 3 a 3 + A§ 

where a 3 = -(ai + a 2 ),b 3 = -(b l + b 2 ),A j = Xj-gBjJ = 1,2,3, B 3 = ~(B 1 + B 2 ), 
rij = 0, 1,2,while A j = ±(lj + 1) are the eigenvalues of euclidean Dirac operator Vq 
on unit sphere with lj = 0,1,2, .... Besides 

Q 1 = \/ A i ~ g 2 a\ ,OL2 = \J A| ,« 3 = \/a§ - 5 2 (ai + a 2 ) 2 • ( 7 . 20 ) 


Further, radial parts of ( 6 . 29 ) arc given at n : j = 0 by ( 7 . 13 ) while at n 3 > 0 by ( 7 . 14 ). 
We have ry = 2 / 3 j-r, where ( 3 j = + p 2 6 2 at j = 1 , 2,3 with & 3 = — (6i + 62), 

Pj = gbj + ( 3 j, j = 1,2, P 3 = -5(61 + 62) + /? 3 , Qj = g 0 + w,, Y) = (ajPj - 
gajujj + gajbj)Qj + g 2 a j b j P j , j = 1,2 , Y 3 = [a 3 /? 3 + c?(ai + a 2 )w 3 - 90:3(61 + 6 2 )]Q3 + 
p 2 (ai + a 2 )(bi + b 2 )P 3 , Z 3 = [(Aj - gBj)Pj + ga 3 go\Q 3 + g 2 a j b j P j , j = 1,2, Z 3 = 
[(A3 + g(-Bi + B 2 ))F , 3 - g(ai + a 2 )go\Q3 + 9 2 (ai + a 2 )(&i + b 2 )P 3 , quantum numbers 
rij arc defined by the relations 

gbiZi — (3\Yi gb 2 Z 2 - (3 2 Y 2 g(b\ +b 2 )Z 3 + fl 3 Y 3 ^ ^ 

" 1 = P1P1Q1 •" 2= ftPjQj '" 3 = --' (7 ' 21) 

Further, Cj of (7.13)-(7.14) is determined from the normalization condition 

roo 1 

(\ Fjl \ 2 + \ Fj2 \ 2 )dr = - . (7.21') 

As a consequence, we shall gain that in (6.29) Vbj e (M 3 ) at any tel and, accordingly, 
T = (\Bi, F 2 . \& 3 ) may describe relativistic bound states with the energy spectrum (7.17)- 
(7.19). 

Finally, it should be noted that the influence of the Dirac-like monopole configura¬ 
tions for gluonic SU(3)-field when Kj / 0 in (5.26) can be treated by the same man¬ 
ner as in Subsection 6.4 if taking cr 2 <frj ~ sin i)4> ? for the eigenspinor <3?j of the twisted 
euclidean Dirac operator P/, on the unit sphere S 2 with the conforming Chern numbers 

k = k\,k 2 , — (ki + k 2 ). 


7.4 iV = 4 

In line with previous Subsection it is not difficult to write out the results for the given case. 
The spectrum is given by uj = uj\ + u> 2 + ui 3 + with 

= u j( n ji (p Aj) = 


-A jg 2 cijbj ± (rij + aj)yj (n 2 + 2 rijaj + A 2 )/qj + g 2 b 2 (n 2 + 2 rijOtj) 

n 2 + 2 rijaj + A 2 - 


1 j — 1, 2,3 , 

(7.22) 
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id 4 — 6 ^ 4 ( 77 . 4 , ( 4 , A 4 ) — 


-A 45 2 a 4&4 ± (n 4 + « 4 )\ J (ra| + 27740:4 + A|)^ + g 2 b \{ n \ + 277404) 


n 


4 + 27740:4 T- A? 


(7.23) 


where a 4 = -(ai + a 2 + a 3 ), 64 = -(61 + 62 + 63 ). A? = Aj - gBj ,j = 1,2,3,4 , B 4 = 
— (B\ + B 2 + /i;!), rij = 0,1,2,.... while Ay = ±((. ; + 1) arc the eigenvalues of euclidean 
Dirac operator Vq on unit sphere with l :j = 0,1, 2,.... Besides 


°tj = yjtf ~ g 2 a 2 j ,j = 1,2,3 ,04 = \J- g 2 (a 1 + a 2 + a 3 ) 2 . (7.24) 

Further, radial parts of (6.29) arc given at rij = 0 by (7.13) while at rij > 0 by (7.14). We 
have rj = 2/3jr, where (3j = \jg g — + g 2 b 2 at j = 1 , 2 , 3,4 with 64 = — (61 + b 2 + b 3 ), 

Pj = Qbj+Pj, j = lj 2 , 3, -P 4 = — <7(&i+£>2+&3)+/24> Qj = Yj = {(XjPj—gajUJj + 

ga.jbj)Qj+g 2 cijbjPj, j = 1,2,3 ,14 = [a4/?4+g(ai+a 2 +a3)w4-£a4(&i+&2+fr3)]Q4 + 
g 2 (ai-\-a 2 + a 3 )(bi+b 2 + a 3 )P 4 , Zj = [(Xj — gBj)Pj + gajfio]Qj + g 2 ajbjPj, j = 1,2,3, 
Z 4 = [(\ 4 -\- g(B\+ B 2 +Bz))P 4 — g(a\ + a 2 + 03 )^ 0 ] Qi+ 9 2 { a i J r a 2 J ro J 3 ){bi+b 2 +b 3 )Pi, 
quantum numbers rij arc defined by the relations 


_ gbjZj — PjYj . _ n 0 0 _ ^(hi + &2 + b%)Z 4 + /34T4 

ni “ PjPjQj ,J ~ ’ ’ ,n4_ ^ 


(7.25) 


Further, Cj of (7.13)-(7.14) is determined from the normalization condition 

fOO 1 

J q (|F jl | 2 + |F j2 | 2 )dr = -. (7.26) 

Thus, in (6.29) 'Tj 6 L^R 3 ) at any t£l and, consequently, T 1 = (4c, ^ 2 , ^3 , T 4 ) may 
describe relativistic bound states with the energy spectrum (7.22)-(7.23). 

Finally, it should be noted that the influence of the Dirac-like monopole configura¬ 
tions for gluonic SU(4)-field when K :I / 0 in (5.28) can be treated by the same man¬ 
ner as in Subsection 6.5 if taking cr 2 &j ~ sin (Ah, for the eigenspinor <1> :I of the twisted 
euclidean Dirac operator P/, : on the unit sphere Er with the conforming Chern numbers 
k = h,k 2 , k 3 , -(ki + k 2 + k 3 ). 


7.5 Remark about the case N > 4 

As is not complicated to see, one can consider general case IV > 4. To obtain the cor¬ 
responding confining solutions for given N one should take a concrete realization of the 
conforming SU(A r )-Lie algebra whose form depends on N, as we have seen above. Then 
the explicit form of the solutions under discussion will also depend on the mentioned re¬ 
alization but it is clear that they will be characterized by real constants aj,bj,Aj,Bj, 
j = 1, •••, N with ajv = — (01 + ci 2 + • • • + njv— 1 )> = — (61 + b 2 + • • • + bjy— 1 ), 

Bn = —{B\ + B 2 + • • • + B n _ 1 ) while one may put A, } = 0. 
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The shape of spectrum for relativistic bound states will obviously have the same form 
as in (7.22) with A j = A j — gBj, rij = 0.1, 2,..., but j = 1,..., N. The corresponding 
modifications for parameters ay, (3j , Pj , Qj of the wave functions in the form (7.13)-(7.14) 
are also evident. The normalization condition (7.26) will contain 1/N in its right-hand 
side. So = (T ], T 2 , • • •, 4 ; ,v) may describe relativistic bound states with the the energy 
spectrum lu = ui\ + u; 2 Y • • • Y ujn- At last, the influence of the Dirac-like monopole 
configurations for gluonic SLJ(A r )-field can be treated by the same manner as in previous 
Sections. 


7.6 Nonrelativistic Limit 

It makes sense to obtain the nonrelativistic limit (i.e. when c —> 00) for spectrum (7.8) 
in order to us to have a possibility of estimating the contribution of relativistic effects. 
Changing a —► a/(he), B —► B/{Tic) and expanding (7.8) in z = 1/c, we get 


iu{n, l, A) = 

±Po c 2 
where 

f{n, A) = 10n 9 A + 120n 7 A 3 + 252n 5 A 5 + 120?r 3 A 7 + 10nA 9 + 

y- ( 'n 10 + 45n 8 A 2 + 210n 6 A 4 + 210n 4 A 6 + 45n 2 A 8 + A 10 ) . (7.28) 

As is seen from (7.27), at c —> 00 the contribution of linear magnetic field to the spec¬ 
trum really vanishes and spectrum becomes in essence purely Coulomb one (modulo the 
rest energy) which corresponds to the lines discussed in Subsection 5.1. 


it 


g 2 a 2 


2h 2 {n T |A|) 2 


A g 2 ab 


rA R n ^ f (rj \ ^ 


z+0{z 2 ), (7.27) 


8 Uniqueness of the Confining Solutions 

8.1 Uniqueness 

Let us consider the question of uniqueness of the confining solutions. The latter ones were 
defined in Section 1 and in Subsection 5.1 as the spherically symmetric solutions of the 
Yang-Mills equations (B3) (at J = 0) containing only the components of the SU(A r )-(ield 
which are Coulomb-like or linear in r. Additionally we impose the Lorentz condition (BA) 
on the sought solutions. As was remarked in Appendix B, the latter condition is necessary 
for quantizing the gauge fields consistently within the framework of perturbation theory 
(see, e. g. Ref. 10), so we should impose the given condition. Let us for definiteness 
consider the case of group SU(3). Under this situation we can take the general ansatz of 
form 


A = p l rdt T A r dr + A^dd + r u Ad(p , 


(8.1) 
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where A$ = r p T and matrices F = a a \ a , A = f3 a \ a , T = 7“A a with arbitrary real 
constants a“, (3 a , 7 “. It could seem that there is a more general ansatz in the form A = 
r Pa a a \ a dt + A r dr + r pa ^ a \ a dd + r Ua /3 a \ a dp but somewhat more complicated consid¬ 
erations than the ones below (see Appendix E) show that all the same we should have 
Ha = H, v a = v for any a so, within the current Subsection, we at once consider this condi¬ 
tion to be fulfilled to avoid unnecessary complications and also we put p a = p at any a for 
simplicity though it is not obligatory (see Appendix E). 

For the ansatz (8.1) the Lorentz condition (BA) takes the form 

<9 r (?’ 2 sin dg rr A r ) + dp(r 2 smdg®® A$) = 0 

which can be rewritten as 


(^(sinrMtf) + sin dd r (r 2 A r ) = 0 , 


( 8 . 2 ) 


wherefrom it follows cot Or x T + d r (r 2 A r ) = 0 while the latter entails 


Ay 


C cot 0r p l T 

r 2 p + 1 


(8.3) 


with a constant matrix C belonging to SU(3)-Lie algebra. Then we can see that it should 
put C = T = 0 or else A r will not be spherically symmetric and the confining one where 
only the powers of r equal to ±1 are admissible. As a result we come to the conclusion 
that one should put A r — A,; = 0 in (8.1). After this we have F = dA + gA A A = 
—p,r p ~ l Tdt A dr + vr v ~^Adr A dip + gr p,+u [r, A]dt A dtp which entails (with the help of 
(A6)) 


*F = pr p+1 sin lOrdi!) A dip - 7 ^~dt Add -— [T, A]dr A dd 

sin d sin d 

and the Yang-Mills equations (B. 3) (at J = 0) turn into 

p(p + 1 )r p+l sin 2 dr = p 2 r p+2l/ [A, [r, A]\ , 

v(v - 1 )r v ~ 2 A = g 2 r 2p,+iy [r, [r, A]] . (8.4) 

It is now not complicated to enumerate possibilities for obtaining the confining solutions in 
accordance with (8.4), where we should put p = — 1, v = 1. 

1 . r= 0 or A= 0. This situation does obviously not correspond to a confining solution 

2. r = CA with some constant C. This case conforms to that all the parameters a" 
describing electric colour Coulomb field and the ones (3 a for linear magnetic colour 
field arc proportional - the situation is not quite clear from physical point of view 
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3. Matrices r, A arc not equal to zero simultaneously and both matrices belong to Car- 
tan subalgebra of SU(3)-Lie algebra. The parameters a a , (3 a of electric and magnetic 
colour fields are not connected and arbitrary, i. e. they should be chosen from exper¬ 
imental data. The given situation is the most adequate to the physics in question and 
the corresponding confining solution is in essence the same which has been obtained 
in Section 2 in the form (2.13)—(2.14). 

One can somewhat generalize the starting ansatz (8.1) taking it in the form A = (r^T + 
A')dt. + (r u A + B')dtp with matrices A' = A a \ a , B' = B a X a and constants A a ,B a . 
Then considerations along the same above lines draw the conclusion that the nontrivial 
confining solution is described by r, A,A',B' belonging to Cartan subalgebra which we 
could already see in the solutions (2.13)—(2.14). 

Clearly, the obtained results may be extended over all SU(A')-groups and even over all 
semisimple compact Lie groups since for them the corresponding Lie algebras possess just 
the only Cartan subalgebra. Also we can talk about the compact non-semisimple groups, 
for example, U(iV). In the latter case additionally to Cartan subalgebra we have centrum 
consisting from the matrices of the form aljy with arbitrary constant a. Really we have 
obtained the confining solutions for U(l)-group in Subsection 2.2. The most relevant phys¬ 
ical cases are of course U(l)- and SU(3)-ones (QED and QCD), therefore we shall not 
consider further generalizations of the results obtained. It should also be noted that the 
nontrivial confining solutions obtained exist at any gauge coupling constant g, i. e. they 
are essentially nonperturbative ones. At last, as we have seen in Sections 5, 6, there exist 
the nontrivial confining solutions containing the Dirac-like monopole parts. The latter are, 
however, not the spherically symmetric ones and we shall here not dwell upon uniqueness 
for such solutions. 

8.1.1 Remark Concerning the Wilson Confinement Criterion 

In our recent paper ED it was shown that the confining solutions under consideration satisfy 
the so-called Wilson confinement criterion formulated as far back as in Ref. f32l (see also 
Ref. Cl). During a long time up to now this criterion remains the dominant one when 
building one or another approach to the confinement problem. As far as is known to us, 
however, so far no explicit solutions of the SU(./V)-Yang-Mills equations obeying the given 
criterion have been found. So, taking into account uniqueness of the confininig solutions 
under discussion, perhaps there exist no other solutions meeting the mentioned criterion. 

8.2 Nonrelativistic Confining Potentials 

As has been mentioned in Subsection 4.1, in meson spectroscopy (see, e.g., Refs. M 
El and references therein) one often uses the nonrelativistic confining potentials. Those 
confining potentials between quarks here are usually modelled in the form a/r + br with 
some constants a and b. It is clear, however, that from the QCD point of view the interaction 
between quarks should be described by the whole SU(3)-field A jL = A“T a , genuinely 
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relativistic object, the nonrelativistic potential being only some component of ,4" surviving 
in the nonrelativistic limit when the light velocity c —> oo. Let us explore whether such 
potentials may be the solutions of the Maxwell or SU(3)-Yang-Mills equations. Though 
this can be easily derived from the results obtained in Subsection 8.1 let us consider the 
given situation directly in view of its physical importance. 


8.2.1 Maxwell Equations 

In the case of Maxwell equations (B. 5) (at J = 0) the ansatz A = A t dt = ( a/r + br)dt 
yields F = dA = (a/r 2 — b)dt A dr, *F = sin d(br 2 — a)dd A dtp and the relation 
d* F = 2 br sin Odr A dd A dtp = 0 entails 6 = 0. 


8.2.2 SU(3)-Yang-Mills Equations 

We use the ansatz 

A = A/X a dt = (A'/r + B'r)dt (8.5) 

with some constant matrices A' = a a X a , B' = /3 a X a ■ Then A/\A = 0 ,F = dA+gAAA = 
dA = (A! /r 2 — B')dtf\dr, *F = sin d {B'r 2 — A')dd A dip, d*F = 2B'r svaddr Add Adp, 
*F A A — /I A *F = —2 [A!, B']rdt A dd A dp. Under the circumstances the Yang-Mills 
equations (B. 3) (at J = 0) are tantamount to the conditions d*F = 0, *FAA—AA*F = 0. 
The former entails B' = 0, then the latter is fulfilled at any A! and we can see that the 
Coulomb-like field A = (. A'/r)dt is a solution of the Yang-Mills equations (B. 3) (at J = 0) 
with arbitrary constant matrix A' which actually has been obtained in Subsection 8.1. In 
principle the ansatz (8.5) might be a solution of (B. 3) with the source of the form 


J = 2 B'r sin ddr A did A dp + 2 g[A!, B']rdt Add A dp = 

*j = *(j“A a dx») = * (—dt + g 9 =^dr) , (8.6) 

' \ r r sm v ) 

but div(j) / 0 and this is not consistent with the only source (5.6) derived from the QCD 
lagrangian (5.1). We can avoid this difficulty putting matrices A !, B' arc not equal to zero 
simultaneously and both matrices belong to Cartan subalgebra of SU(3)-Lie algebra. Then 
B' = 4sA 3 + /4iAs and for consistency with the only admissible source of (5.6) we should 
require source of (5.6) to be equal to one of ( 8 . 6 ) which entails 


g^( ltl ^I 3 )X a ^X a dx^ 


2((3 3 X 3 + P 8 X 8 ) dt 

r 


wherefrom one can conclude that 

_ _ 2 / 3 3 

g^^t ® h)^ a ^ = 0, a + 3 ,8 , g^^t ® = — , 

r 


_ 20 s _ 

9^{it <8>/ 3 )A 8 4' = 0- , 5 ^( 7 ^ ®/ 3 )A a 'I' = 0,a = 1, ..., 8 , /J ^ t , 


(8.7) 
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which can obviously be satisfied only at /? 3 r '-' (3^ r '~ / 'T —> 0 at each point of Minkowski 
spacetime, i. e., really matrix B' = 0 again. All the above can easily be generalized to any 
N > 1. 

As a result, the potentials employed in nonrelativistic approaches do not obey the 
Maxwell or Yang-Mills equations. The latter ones are essentially relativistic and, as we 
have seen, the components linear in r of the whole A fl are different from A t and related 
with magnetic (colour) field vanishing in the nonrelativistic limit. 

8.3 Remark about Search for Nonrelativistic Confining Potentials 

The above results make us cast a new glance at search of many years for nonrelativistic 
potentials modelling the confinement. Many efforts were devoted to the latter topic, for 
example, within the framework of lattice gauge theories or potential approach (see, e.g., 
Refs. El and references therein). It should be noted, however, that almost in all literature 
on this direction one does not bring up a question: whether such potentials could (or should) 
satisfy the Yang-Mills equations? As is clear from the above the answer is negative. That is 
why the mentioned approaches seem to be inconsistent. 


9 Application to the Charmonium Spectrum 

As we have emphasized in Subsection 5.1, all considerations in the given paper develop 
with the aim of further physical applications or else the results obtained could be only of 
academic interest. Though there are no obstacles to apply the results (in the most phys¬ 
ically interesting case N = 3) to any meson, up to the moment of writing this paper all 
applications were concentrated on quarkonia (charmonium and bottomonium) and the re¬ 
sults are contained in Refs. mmm. Referring for more details to those papers, we shall 
here outline only the main conclusions drawn from the mentioned papers which confirm 
the physical picture underlying considerations of the given paper and warrant the linear 
confinement scenario described in Subsection 4.2. 

9.1 Relativistic Spectrum of Charmonium 

We can use (7.17)—(7.19) with various combinations of signes (±) before second summand 
in numerators of those formulas. In Refs. mmm due to some reasons (inessential now) 
the combination (+ 4—) was employed and besides the replacement U 2 —■ y w';j was made. 
Let us rewrite, e. g., the results out of Ref. 1351 about spectrum according to (7.17)—(7.19) 
but keeping the mentioned combination (+ 4—). Then numerical results for constants 
parametrizing the charmonium spectrum are shown in Table 1. 

One can note that the obtained mass parameter /io is consistent with the present-day 
experimental limits ESI where the current mass of c-quark (2//q) is accepted between 1.1 
GeV and 1.4 GeV. As for the gauge coupling constant g then its value has been chosen 
in accordance with many recent considerations 1371 . wherefrom one can conclude that the 
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Table 1: Gauge coupling constant, mass parameter //(> and parameters of the confining 
SU(3)-gluonic field for charmonium. 


9 

go 

(GeV) 

ai 

a 2 

bi 

(GeV) 

62 

(GeV) 

Bi 

b 2 

0.46900 

0.62500 

2.21104 

-0.751317 

20.2395 

-12.6317 

6.89659 

6.89659 


strong coupling constant a s = g 2 is of order 0.22 ~ 0.469 2 at the scale of the c-quark 
current mass. At last, as to parameters of solution (2.13), then they only shift the 
origin of count for the corresponding energies and we can consider A\ = A 2 = 0, as was 
mentioned in Subsections 6.2, 6.4. 

With the constants of Table 1 the present-day levels of charmonium spectrum were cal¬ 
culated with the help of (7.17)—(7.19) so Table 2 contains experimental values of these levels 
(from Ref. l36ll ) and our theoretical ones computed according to the shown combinations 
(we use the notations of levels from Ref. (361). 


Table 2: Theoretical and experimental charmonium levels. 


State 

3 

Theoret. energy ej = ^ uok 

k=i 

(GeV) 

Experim. 

value 

(GeV) 

fic(l S) 

£i = wi(0,0, -1) + w 2 ( 0 ,0, -1) + w 3 (0,0, -1) = 2.979597 

2.979600 


e 2 = wi(0,0, -1) + w 2 (0,0, -1) + w 3 (0,0,1) = 3.096913 

3.096916 

Xco(l P) 

£3 =wi(0,0,-1)+ w 2 (0,0,1) + w 3 (0,0,-1) = 3.415186 

3.415190 

Xci(l P) 

£4 = wi(0,0,1) + w 2 (0,1, -1) + w 3 (2, 0,1) = 3.505304 

3.510590 

hc(lP) 

£5 = wi(0,0, -1) + wa(0,0,1) + w 3 (0,0,1) = 3.532503 

3.526210 

Xc 2 (l P) 

£ 6 =wi(0,1,-1) + w 2 (l, 1,-1) + w 3 (l, 1,-1) = 3.553097 

3.556260 

Vc(2S) 

£7 = wi(0,0,1) + wa(0,1, -1) + ws(l, 0,-1) = 3.671608 

3.65400 

V>(2 S) 

£ S = wi(0,1, -1) + W 2 (l, 1, -1) + w 3 (2, 1,1) = 3.674025 

3.685093 

^(3770) 

£9 = wi(0,0,1) + W2(0,0,1) + w 3 (2, 0, -1) = 3.775598 

3.770000 

V>(3836) 

£10 =wi(0,1,-1)+ u 2 (0,1,1) +w 3 (0,0,1) = 3.833640 

3.836000 

X(3872) 

£n = wi(0,1,-1) + w 2 ( 0,1,1) + w 3 (0,1,1) = 3.871672 

3.872000 

^(4040) 

£12 = wi(0,0,1) + w 2 (0,1,-1) + W 3 (l, 1,1) = 4.042660 

4.040000 

V’(4160) 

£i 3 =wi(0,0,-1)+ wa(0,1,1)+ w 3 (0,0,-1) = 4.153765 

4.159000 

i/>(4415) 

£14 =wi(0, 0,-1)+ w 2 (l, 0,-1)+ w 3 (2, 1,1) = 4.409260 

4.415000 
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9.2 Nonrelativistic Limit 

One can be interested in estimating the contribution of relativistic effects, i. e. those con¬ 
nected with magnetic colour field linear in r. This has been done in Refs. 00 with the 
help of relations (7.27)-(7.28). The contribution of relativistic effects can amount to tens 
per cent and they cannot be considered small. Moreover, the more excited the state of char- 
monium the worse the nonrelativistic approximation. The physical reason of it is quite clear. 
Really, we have seen in the nonrelativistic limit (see the relations [7.27)-(7.28)] that the pa¬ 
rameters 612, -61,2 of the lineal - interaction between quarks vanish under this limit and the 
nonrelativistic spectrum is independent of them and is practically getting the pure Coulomb 
one. As a consequence, the picture of linear confinement for quarks should be considered 
an essentially relativistic one while the nonrelativistic limit is very crude approximation. 
In fact, as follows from exact solutions of SU(3)-Yang-Mills equations of (2.13), the lin¬ 
ear interaction between quarks is connected with magnetic colour field that dies out in the 
nonrelativistic limit, i.e. for static quarks. Only for the moving rapidly enough quarks the 
above field will appear and generate linear confinement between them. So the spectrum 
will depend on both the static Coulomb electric colour field and the dynamical magnetic 
colour field responsible for the linear confinement for quarks which is just confirmed by the 
relations (7.17)—(7.19). In our case, the interaction effect with the magnetic colour field is 
taken into consideration from the very outset which just reflects the linear confinement at 
large distances. 

9.3 Electromagnetic Transitions 

We can specify the obtained above charmonium spectrum. The fact is that the relations 
(7.17)—(7.19) permit various parametrizations of the charmonium spectrum (see Refs. 0 
0) and therefore it should impose further conditions to fix a certain parametrization among 
several possible ones. For example, one can compute widths of elecromagnetic transitions 
among levels of charmonium, in particular for transitions J/'ip(\S) —>■ r] c (lS) + 7 and 
X c o(l-P) — * J/ip{lS) + 7. In Ref. l35l the widths of the latter transitions have been cal¬ 
culated in dipole approximation that allowed one to use the corresponding wave functions 
described in Subsection 7.3. Results of computation supplies us with an additional justifi¬ 
cation for the choice of parameters of the SU(3)-confining gluon field adduced in Table 1 
and allows us to conclude that dipole approximation is not enough for the second transition 
of the ones under discussion. The question now is what gluon concentrations are in the 
mentioned SU(3)-confining gluon field. 

9.4 Estimates of Gluon Concentrations and Magnetic Colour Field Sstrength 

To obtain necessary estimates we shall use an analogy with classical electrodynamics where 
is well known (see e. g. 03) that the notion of classical electromagnetic field (a photon 
condensate) generated by a charged particle is applicable only at distances > > the Comp¬ 
ton wavelength A c = 1/m for the given particle. If denoting A b the de Broglie wavelength 
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of the particle then A b = 1 /p with the relativistic impulse p = rnv /\/l — v 2 while v 
is the velocity of the particle (as a result, A c = Ar at v = 1/ \/2) so one can rewrite 
A b = A c Vl — v 2 /v and it is clear that A# —> 0 when v —» 1 (ultrarelativistic case), i. e., 
the particle becomes more and more point-like one. Accordingly, one can conclude that in 
the latter case the notion of classical electromagnetic field generated by a charged ultrarela¬ 
tivistic particle is applicable at distances >> A b- Under the circumstances, if the ultrarela¬ 
tivistic charged particle accomplishes its motion within the region with characteristic size of 
order tq then in the given region the electromagnetic field generated by the particle may be 
considered as classical one at ro >> A b- For example, in the case of positronium we have 
ro ~ 2ao, where ao ~ 5.29 • 10 -11 m is the Bohr radius, so ro >> A e ~ 3.86 • 10 _1! nr, 
the electron Compton wavelength, i. e., the electric Coulomb interaction between electron 
and positron in positronium can be considered classical electromagnetic field. 

Passing on to QCD, gluons and quarkonia, it should be noted that quarks in quarkonia 
accomplish a finite motion within a region of order 1 fm = 10“ 15 m. Then, as is seen from 
the radial parts of the wave functions described in Subsection 7.3, the quantity l//3j permits 
to be considered a characteristic size of the j-th colour component of the given quarkonium 
state and, consequently, we can take the magnitude 

(9 - 1) 

for a characteristic size of the whole quarkonium state and, in line with the above, we should 
consider the confining SU(3)-gluonic Yang-Mills field of (2.13) or (2.14) to be classical 
one when ro >> Ar, the de Broglie wavelength of the corresponding quarks forming 
quarkonium. 

On the other hand, a classical electromagnetic field (photon condensate) conforms to 
the large photon concentrations for every frequency presented in the field 1251 . Then in 
QCD we should require the large gluon concentrations in the given classical gluonic field 
(gluon condensate). For the necessary estimates we shall employ the Too-component of the 
energy-momentum tensor for a S U(A r )-Yang- M i 11 s field 



To estimate the given concentrations we can employ Too-component (volumetric energy 
density) of the energy-momentum tensor of (9.1') and, taking the quantity u = T, the 
whole decay width of the quarkonium state, for the characteristic frequency we obtain the 
sought characteristic concentration n in the form 


n = 


2qq 

r 


(9.2) 


It is not complicated to obtain the curvature matrix (field strentgh) corresponding to the 
solution (2.13) or (2.14) 


F = T“ A a TzT A dx u = — 


d r (A^\ a )dt A dr + d r {A^X a )dr A dp 


(9.3) 
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which entails the only nonzero components 

fli — «2 j-,8 (ai + a2)\/3 


El3 _ “1 “2 ^8 _ TU^VO ,-,3 

± f r — ~ ~ i r+*. — ~ ~ i -t — 


2 r 2 


2r 2 ’ tr 

and, in its turn, 

Too = T tt = ^ [(T 3 ) 2 + (Ff r ) 2 


b\ — &2 p8 _ (bi + b 2 )VS 

5 -*■ ry? 2 ? " Ty> 2 


(9.4) 


+ 


1 


4r 2 sin 2 $ L rLp 


«) 2 + ) 


f + Or a.2 + cl 2 b 2 + 6162 4 - ^2 


+ 


167T ^ r 4 ' 3r 2 sin 2 d y ’ 

so, further putting sin 2 d = 1/3 for simplicity, we can rewrite (9.5) in the form 


(9.5) 


rri _ rr-iCOlll 1 /Tnlill 

1 00 — 1 00 + 1 00 


(9.6) 


conforming to the contributions from the Coulomb and linear parts of the solutions (2.13) 
or (2.14). The latter gives the corresponding split of n from (9.2) 


n — n cou i + nii n . 


(9.7) 


Using the Hodge star operator in 3-dimensional euclidean space [where ds 2 = g llv dxJ l 0 
dx v = dr 2 + r 2 (dd 2 +sin 2 dd(p 2 ), *(dr Add) = sin d dip,* (dr A dip) = —dd/sind,*(ddA 
dip) = dr/(r 2 sin 0 1, we can confront components F)L. F/ p with 3-dimensional 1-forms of 
the magnetic colour field 


H 3 = 


bi - b 2 
2sinr? 


dd, H 8 


(bi + b 2 )V3 M 

2 sin d 


that are modulo equal to \J g^H 3,8 H 3 ’ 8 or 

h 3 = \bi ~ b 2 \ H s = \bi + b 2 \V3 
2rsinr9 ’ 2rsin$ 


(9.8) 


where we also consider sin 2 d = 1/3 for simplicity. 

For comparison we shall also estimate the photon concentration in the ground state of 
the positronium. As is known historically j38j, the analogy between the latter system and 
quarkonia played the important part when building the quarkonia models. For positronium 
we have the electromagnetic Coulomb interaction A = A t dt = (e/r)dt which entails 
F = Ftrdt A dr = ( e/r 2 )dt A dr and 


Too = T tt = 



Oiem 

47T?’ 4 


(9.9) 


with a em = e 2 = 1/137.0359895. 
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9.5 Numerical Results and Concluding Remarks 

When computing for the ground state of charmonium we used its present-day whole decay 
width r = 17.3 MeV 1361 . while the calculation ro of (9.1) gives tq = r\ (see Table 3). 

In the positronium case we employed the widths i~o = 1/ro (parapositronium) and 
Ci = 1/n (orthopositronium), respectively, with the life times tq = 1.252 • 10 -10 s, 
t i = 1.377 • 1CT 7 s 1.381 while ro = 2ao with the Bohr radius ao = 0.529177249 • 10 5 fm 

ESI. 

Tables 3, 4 contain the numerical results for both the cases, where, when calculating, 
we applied the relations 1 GeV -1 « 0.21030893 fm , 1 T ~ 0.692508 x 10 -15 GeV 2 . 

Table 3: Gluon concentrations and magnetic colour field strengths in the ground state of 
charmonium. 


ric(isy. 

n = 0.0399766 fm 








r 

(fm) 

n coul 

(m“ 3 ) 


^lin 

(m —3 ) 


n 

(m —3 ) 


H 3 

(T) 


ff 8 

(T) 


O.lri 

0.550649 

10 58 

0.727630 

1 0 55 

0.551377 

1 0 58 

0.216259 

10 19 

0.866919 

10 18 

r l 

0.550649 

10 54 

0.727630 

10 53 

0.623412 

10 54 

0.216259 

10 18 

0.866919 

10 17 

10ri 

0.550649 

10 5 ° 

0.727630 

10 51 

0.782695 

10 51 

0.216259 

10 17 

0.866919 

10 16 

1.0 

0.140637 

Ol 

O 

0.116285 

10 51 

0.117691 

10 51 

0.864530 

10 16 

0.346565 

10 16 

a 0 

0.179347 

o 

CO 

o 

0.415260 

10 41 

0.415260 

10 41 

0.163372 

10 12 

0.654912 

10 11 


Table 4: Photon concentrations in the ground state of positronium. 


r 0 = 2a 0 = 2 ■ 0.529177249 ■ 10 s fm 


Parapositronium 

Orthopositronium 

r 

TTpara 

T^ortho 

(fm) 

(m- 3 ) 

fm" 3 ) 

O.Olro 

0.888025 ■ 10 46 

0.976685 ■ 10 49 

0.1r 0 

0.888025 ■ 10 42 

0.976685 ■ 10 45 

’’O 

0.888025 ■ 10 38 

0.976685 ■ 10 41 

2r 0 

0.555015 ■ 10 37 

0.610428 ■ 10 4 ° 


Then, as is seen from Tables 3, 4, qualitative behaviour of both the concentrations is 
similar. At the characteristic scales of each system the concentrations are large and the 
corresponding fields (electric and magnetic colour ones or electric Coulomb one) can be 
considered classical ones. For charmonium the part ro cou i of gluon concentration n con¬ 
nected with the Coulomb electric colour field is decreasing faster than nu n , the part of n 
related to the linear magnetic colour field, and at large distances ny m becomes dominant. 
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Under the circumstances, as has been said above, we can estimate the quark velocities in 
the charmonium state under discussion from the condition 



with the c-quark Compton wavelength \ q = l/(2//-o) ~ 0.168247 fm so, taking the de 
Broglie wavelength A b = O.lri with r\ from Table 3, we obtain v ~ 0.999718, i. e., the 
quarks in charmonium should be considered the ultrarelativistic point-like particles. This 
additionally confirms the conclusion of Refs. GEl that the relativistic effects are extremely 
important for the confinement mechanism. As a result, the confinement scenario described 
in Subsection 4.2 may really occur. At last, we can see that strength of magnetic colour 
field responsible for linear confinement reaches huge values of order 10 1 '-10 18 T. For 
comparison one should notice that the most strong magnetic fields known at present have 
been discovered in magnetic neutron stars, pulsars (see, e. g., Ref. (39)) where the corre¬ 
sponding strengths can be of order 10 9 -10 10 T. So the characterictic feature of confinement 
is really the very strong magnetic colour field between quarks which we have emphasized 
in Subsection 5.1. In a certain sense the essence of confinement can be said to be just in 
enormous gluon concentrations and magnetic colour field strentghs in space around quarks. 


10 Conclusion 

Throughout the paper we moved step-by-step forward in analysing the Dirac-Yang-Mills 
system of equations (5.4)—(5.5) derived from the QCD lagrangian. The aim we pursued was 
to obtain a scenario for linear confinement of quarks, at any rate, in mesons and quarkonia. 
It seems to us we succeeded in finding the suitable quantitative description for the given 
phenomenon. As we could see, crucial step here consisted in studying exact solutions of the 
SU(3)-Yang-Mills equations modelling confinement and the corresponding modulo square 
integrable solutions of the Dirac equation. Techniques of finding the mentioned solutions 
were based to a great degree on those borrowed from black hole theory. In this respect 
our approach is conventional for physics whose whole history shows that very often the 
methods developed for solving some problems proved to be extremely useful for analysing 
a number of other ones in a perfectly different region of physics. 

At the end of our considerations let us summarize the main features of the confinement 
mechanism developed in the given paper. 

1. The whole approach is based on the exact solutions of the Yang-Mills equations and 
on the corresponding modulo square integrable solutions of Dirac equation. The so¬ 
lutions in question are essentially unique ones: for the confining solutions of the 
SU(3)-Yang-Mills equations that was shown in Section 8 whereas uniqueness for the 
conforming solutions of Dirac equation (5.4) follows from the ansatz 7.1'. Namely, 
the mentioned ansatz leads to (7.3) which is the confluent hypergeometric equation 
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and (as is known from literature on special functions, see, e. g., Ref. (H) the latter 
just possesses the only suitable solutions for us in the form of Laguerre polynomials 
which in essence determines the spectrum of relativistic bound states [see relations 
(7.7)—(7.8)]. Consequently, we found the only compatible solutions of the system 
(5.4)—(5.5) which can have pretensions of describing the confinement mechanism. 
Another matter that lagrangian QCD probably allows one to develop some other ap¬ 
proaches [not based on compatible solutions of (5.4)—(5.5) ] to confinement but our 
one seems to be the most natural since it is practically identical to the standard ap¬ 
proach of quantum mechanics and QED to description of bound states in external 
electromagnetic fields. In other words, this approach should have been one of the 
very first approaches as soon as the QCD lagrangian was written. Historically, how¬ 
ever, this way was rejected due to incomprehensible reasons. 

2. Two main physical reasons for linear confinement in the mechanism under discus¬ 
sion are the following ones. The first one is that gluon exchange between quarks 
is realized with the propagator different from the photon one and existence of such 
a propagator is direct sequence of the unique confining solutions of the Yang-Mills 
equations. The second reason is that, owing to the structure of mentioned propaga¬ 
tor, gluon condensate (a classical gluon field) between quarks mainly consists of soft 
gluons ( see Subsection 4.1) but, because of that any gluon also emits gluons, the 
corresponding gluon concentrations rapidly become huge and form the linear con¬ 
fining magnetic colour field of enormous strengths which leads to confinement of 
quarks. Under the circumstances physically nonlinearity of the Yang-Mills equations 
effectively vanishes so the latter possess only the unique confining solutions of the 
abelian-like form (with the values in Cartan subalgebra) that describe the gluon con¬ 
densate under consideration. Moreover, since the overwhelming majority of gluons 
are soft they cannot leave hadron (meson) until some gluon obtains additional energy 
(due to an external reason) to rush out. The latter seems to be observable just in 
the so-called 3-jets events (for more details see, e. g., Ref. El). So we deal with 
confinement of gluons as well. 

3. The approach under discussion equips us with the explicit wave functions that is 
practically unreachable in other approaches, for example, within framework of lat¬ 
tice gauge theories or potenial approach. Namely, for each meson there exists its own 
set of real constants aj , Aj , bj, Bj parametrizing the mentioned confining gluon field 
(the gluon condensate) and the corresponding wave functions while the latter ones 
also depend on /j,q, the reduced mass of the current masses of quarks forming me¬ 
son. It is clear that constants aj, Aj, bj, Bj,no should be extracted from experimental 
data. This circumstance gives possibilities for direct physical modelling of internal 
structure for any meson and for checking such relativistic models numerically. 

4. Finally, there is also an interesting possibility of indirect experimental verification 
of the confinement mechanism under discussion. Really solutions (2.6) point out 
the confinement phase could be in electrodynamics as well. Though there exist no 
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elementary charged particles generating a constant magnetic field linear in r, the 
distance from particle, after all, if it could generate this elecromagnetic field config¬ 
uration in laboratory then one might study motion of trial charged particles in that 
field. The confining properties of the mentioned field should be displayed at classical 
level too but the exact behaviour of particles in this field requires certain analysis of 
the corresponding classical equations of motion. 

At this point we would like to mark the end of the given paper on pinnig our hopes on 
studying confinement with subsequent development of the results obtained here, in the first 
place, by further applications to concrete mesons. 

11 Appendix A: Hodge Star Operator * on Minkowski 
Spacetime in Spherical Coordinates 

Let M is a smooth manifold of dimension n so we denote an algebra of smooth func¬ 
tions on M as F(M). In a standard way the spaces of smooth differential p-forms A P (M) 
(0 < p < n) are defined over M as modules over F(M) (see, e. g. Refs. t20V ). If a 
(pseudo)riemannian metric G = ds 2 = g lw dx IJ <g> dx u is given on M in local coordinates 
x = (x 1 ') then G can naturally be continued on spaces A P (M) by relation 

G(a, ft) = det{G(ai, (3j)} (A.l) 

(j) 

fora = a\/\a 2 ...f\a p , [3 = (3\ A/3 2 ... A/? p , where for 1-forms a* = a.y dx^ , f3j = [31 dx u 
we have G(ai,(3j ) = g^ctfi) (3^ with the Cartan’s wedge (external) product A. Under 
the circumstances the Hodge star operator *: A P (M) —► A n ~ p (M) is defined for any 
a € A P (M) by 

a A (*a) = G(a,a)ujg (A.2) 

with the volume n-form u> g = \J\ d(-'X(g /lu )\dx 1 A ...dx n . For example, for 2-forms F = 
Fgudx 11 A dx u we have 

F A *F = (jf a y 0 - g va )F^F a pV5 dx 1 A dx 2 • • • A dx n , y, < is, a < (3 (A. 3) 

with 5 = | det(p /i!y )|. If s is the number of (-1) in a canonical presentation of quadratic 
form G then two most important properties of * for us are 

* 2 = (_i)p(»-p)+« , (A4) 

*{fiai + f 2 a 2 ) = fi(*ai) + f 2 (*a 2 ) (A. 5) 

for any fi, f 2 £ F[M), ai,a 2 € A P (M), i. e., * is a F(M )-linear operator. Due to (A.5) 
for description of ^-action in local coordinates it is enough to specify ^-action on the basis 
elements of A P (M), i. e. on the forms dx 11 A dx 12 A ... A dx lp with i\ < i 2 < ... < i p 
whose number is equal to C(’ = ■ 
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The most important case of M in the given paper is the Minkowski spacetime with 
local coordinates t, r, i ). ip, where r, id, ip stand for spherical coordinates on spatial part of 
M. The metric is given by (1.1) and we shall obtain the ^-action on the basis differential 
forms according to (A.2) 

*dt = r 2 sin iddr A did A dip, *dr = r 2 sin ddt A did A dip, 

*did = — r sin iddt A dr A dip, *dip = rdt A dr A did, 

*{dt A dr) = —r 2 sin iddid A dip ,*(dt A did) = sin iddr A dip , 

*(dt A dip) = - -dr A did , *(dr A did) = sin iddt, A dip , 

sin id 

*(dr A dip) =- -dt A did , *(did A dip) = - -dt A dr , 

sin id r z sin id 

*(dt A dr A did) = -dip , *(dt A dr A dip) =- -did, 

r r sin id 

*(dt A did A dip) = - -dr, *(dr A did A dip) = - -dt, (^4-6) 

r z sin 1 / r z sin id 

so that on 2-forms * 2 = — 1, as should be in accordance with (AA). 

At last it should be noted that all the above is easily over linearity continued on the 
matrix-valued differential forms (see, e. g., Ref. ED),i- e., on the arbitrary linear combina¬ 
tions of forms aj 1 j 2 ...j p (ix n A dx 12 A ... A dx' v , where coefficients belong to some 

space of matrices V, for example, a SU(A r )-Lic algebra. In the latter case, if T a are matri¬ 
ces of generators of the SU(A r )-Lie algebra in iV-dimensional representation, we continue 
the above scalar product G on the SU(iV)-Lie algebra valued 1-forms A = A" T„dx 1 ’ and 
B = B^Tbdx v by the relation 

G(A, B) = g^A;B b u Ti(T a T b ), (A.7) 

where Tr signifies the trace of a matrix, and on linearity with the help of (A.l), (A.7) can 
be continued over any SU(A r )-Lie algebra valued forms. 

12 Appendix B: SU(A)-Lie Algebras and their Cartan 
Subalgebras 

SU(A^)-Yang-Mills Fields as Connections in Vector Bundles 

The most convenient mathematical treatment of classical SU(A r )-Yang-Mil Is fields is, to 
our mind, the one decribing those fields as connections in vector (or, which is equivalent, in 
principal) bundles over some manifold M. When M is the Minkowski space the situation 
is simplified so long as all the bundles over M with topology K 4 are trivial. Referring for 
more details, e. g. to Refs. (40], we shall restrict ourselves to a few remarks. 
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If denoting £ the standard (trivial) iV-dimensional vector bundle over Minkowski space- 
time we can introduce an SU(iV)-connection (a classical SU(iV)-Yang-Mills field) in £ as a 
SU(A r )-Lie algebra valued form (the connection matrix) A = A^dxJd A fl = A“T a while 
the matrices T a form a basis of the Lie algebra of group SULY) in iV-dimensional space 
(we consider T a hermitian, as is accepted in physics), a = 1, N' 2 — 1. Then the curvature 
matrix (field strentgh) for ^-bundle is F = dA + gA A A = F" w T a dx^ A dx u with the exte¬ 
rior differential d (for example, d = dfdt + d r dr + d^d:0 + d v dip in coordinates t, r, i), tp) 
while a gauge coupling constant g is introduced from physical considerations. Any smooth 
function (the gauge transformation) S : M —> SU(iV) gives the so-called trivialization of 
^-bundle and if As and Fs are the connection and curvature matrices, respectively, for 
the given trivialization then they are related to the previous ones by relations (the gauge 
transformations) 

A s = S' 1 AS + -S- l dS, F s = 5 _1 F5, (BA) 

g 

where the factor i is taken in order to an SU(A r )-Lic algebra could be chosen from the 
hermitian matrices. 

Mathematically A and F are linked by the Bianchi identity holding true for any con¬ 
nection 

dF = FAA — AaF, (B. 2) 

and physical considerations impose the Yang-Mills equations 

d * F = g(*F A A — A A *F) + J , (B .3) 

where * means the Hodge star operator conforming to a Minkowski metric, for instance, in 
the form of (1.1), while J is a source, i. e., some differential SU(A r )-algcbra valued 3-form. 
An arbitrary SU(A' r )-connection does not obey the equations (B. 3) except for the so-called 
self-dual fields for those F = *F and the equations (B. 2) and (B. 3) become the same (at 
J = 0, g = 1). But, as was remarked in Appendix A, in Minkowski spacetime * 2 = — 1 on 
2-forms that entails F = 0 for self-dual fields. 

Also it should be noted that the solutions of (B. 3) arc usually believed to obey an 
additional condition. In the present paper we take the Lorentz condition that can be written 
in the form div(A) = 0, where the divergence of the Lie algebra valued 1-form A = 
A^Tadx 11 is defined by the relation (see, e. g. Refs. 12011 ) 

div(A) = -L d,(Vdg^A v ) . (BA) 

It should be emphasized that the Lorentz condition is necessary for quantizing the gauge 
fields consistently within the framework of perturbation theory (see, e. g. Ref. Q), so we 
should impose the given condition. 

At last, the case of group U(l) is also very important. We then have the case of clas¬ 
sical electromagnetic field (the corresponding Lie algebra consists from real numbers), the 
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Bianchi identity is converted into the the first pair of Maxwell equations dF = 0 while the 
Yang-Mills equations {B. 3) become the second pair of Maxwell equations 

d*F = J (B.5) 

with F = dA, A = A^dx 1 ' and J = j fJ * (dx 11 ) with a 4-dimensional electromagnetic 
density current j = j^dx 1 ', where = j/Jx) arc some functions on K 4 . 

Let us describe the explicit realizations of SU(iV)-Lie algebras by hermitian matrices 
that arc needed for our aims and also indicate the so-called Cartan subalgebras in them 
which is important in the main paid of the paper. By definition, a Cartan subalgebra is a 
maximal abelian subalgebra in the corresponding Lie algebra, i. e., the commutator of any 
two matrices of the Cartan subalgebra is equal to zero. Dimension of Cartan subalgebra (as 
a vector space) for SU(A r )-Lic algebra is equal to N — 1. 

N = 2 

In this case we can take T a = cr a at a = 1, 2,3, where a a arc the ordinary Pauli matrices 



We shall also notice that 

Tv(a k aj) = ^5 kj . (B.7) 

The Cartan subalgebra is generated by 03 . 


N = 3 

In the given situation we can take T a = \ a , where X a arc the Gell-Mann matrices 

/0 1 0\ /0 —i 0\ /I 0 0\ 

Ai = 1 0 0 ,A 2 = * 0 0 ,A 3 = 0 -1 0 , 

Vo 0 0/ Vo 0 0/ Vo 0 0/ 

/O 0 1\ /o 0 -*\ /o 0 0\ 

a 4 = 0 0 0 , a 5 = 0 0 0 , a 6 = 0 0 1 , 

Vi 0 0/ \i 0 0 J Vo 1 0 ) 

(B. 8) 


(B. 9) 


/0 0 0 \ /I 0 0 \ 

a 7 = 0 0 -i ,a 8 = -= 0 1 0 , 

Vo i 0 / v 3 Vo o - 2 / 

and 

Tr(A fe Aj) = 2 S kj . 

The Cartan subalgebra is generated by A 3 and As- 
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N = 4 


Following Ref. ATI we take (I 2 stands for the unit matrix 2x2) 


T\ 




,t 2 






t 4 




( 0 —il 2 \ (I 2 0 \ 

U o; - rs - f o -ij 


(B. 10 ) 


and the products of matrices / 7 — I l \. T§ — 7 ) / 3 , Tg — T 4 Tg, T\o — T 2 Tj. / ; — T 2 T§, 
T 12 = T 2 T 6 , T13 = T 3 T 4 , T 14 = T 3 T 5 , T15 = T 3 T 6 for the rest of the Lie algebra basis so 
the Cartan subalgebra is generated by T 3 , Tg, T15 with 


T '» = T = r «=(o -"3) 


(£• 11 ) 


Then 

TV(T fc T,) = 4 <y fcj . ( 5 . 12 ) 


13 Appendix C: Eigenspinors of the Euclidean Dirac Operator 

on § 2 at A = ±1 


When using the wave functions obtained in Section 7 in any applications (see, e.g., Ref. 
{ 35 )) one needs the explicit form for eigenspinors of the euclidean Dirac operator Vq (see 
Subsection 6 . 1 ). Though it is in general given by the relation (6.6), for applications, as a 
rule, it is sufficient to restrict oneself to small eigenvalues of Vq. Let us write out, therefore, 
the eigenspinors of the euclidean Dirac operator Vq corresponding to the eigenvalues A = 
±1 in explicit form. 

If A = ±(l + 1 ) = ±1 then l = 0 and from (6.6) it follows that k = l + 1/2 = 1 / 2 , 
rn '| < 1/2 and we need the functions P±i/ 2 ±i/ 2 - But according to Ref. f 27 l we have 

PL = cos 2fc (tf/ 2 ), Ptk = * 2fc sin 2fc (r?/ 2 ) and, besides, = Pl k _ k , P k l kk = P^_ k 

that entails the eigenspinors for A = — 1 in the form 


C / cos | + i sin | 
2 V cos 4 — i sin f 


e^ /2 ,$ 


C( cosf + zsinf \ i(p/2 

2 V — cos | + i sin f ) 


while for A = 1 the conforming spinors are 


C / cos | — i sin | 
2 v cos | + i sin | 


e iv/ 2 ,<f> 


C / -cos | + 7 sinf \ i(p/2 

2 V cos | + i sin | ) 


(C.l) 


( 67 . 2 ) 


with the coefficient C = a/ 1 /( 27 t). 
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14 Appendix D: Condition <72 ~ sin$<I> 


Let us compute spinor 4>o = (a 2 — sin d) ( I> where, for example, $ is the first spinor in 
((7.1). We obtain 


<Ln = 


d 


d 


— sin — (2 + cos d) I, + i cos — (2 — cos d) 

2 VI/ 2 V 1 


-1 


J<p /2 


(77.1) 


so that both components of <1+ are modulo 


1 + 3 sin^ d 


87T 


< 


0.4 < 1 


i. e. condition o+b ~ sin 0 <I> is good enough fulfilled. This holds true for all the eigen- 
spinors of Vq, or, more generally, of V n which we shall not dwell upon. 


15 Appendix E: More General Ansatz 

As was remarked in Subsection 8.1, there can be a more general ansatz than (8.1) for finding 
the confining solutions, namely it is 

A = r^ a a a \ a dt + A r dr + A,$dd + r Ua /3 a X a dp , (77.1) 


with A$ = r pa j a X a and arbitrary real constants a", (3 a , 7+ Under this situation the Lorentz 
condition (BA) for the given ansatz entails cot 0r p "p°X a + d r (r 2 A r ) = 0 wherefrom it 
follows 


A r 


C_ 

+2 


cot dr pa 1 
Pa + 1 


I a \a 


(77.2) 


with a constant matrix C belonging to SU(3)-Lie algebra. Then we can see that it should 
put C = 7“ = 0 or else A r will not be spherically symmetric and the confining one where 
only the powers of r equal to ±1 are admissible. As a result, we come to the conclusion 
that one should put A r = A$ = 0 in (77.1). 

After this we have 


cL4 = — 1 a a X a dt A dr + v a r Va 1 /3 a X a dr A dp, 

F = dA + gA A A = dA + ga a /3 b r pa+Ub [X a , Xb]dt A dp 

and, with the help of (A.6), we obtain 


*F = p a r Pa+1 a a X a sin ddd A dp 


v a r Ua 1 (3 a X a 
sinr? 


dt A dd— 


gr 


^ +Ub a a f3 b 


[Aaj Aft] 


dr Add 


(77.3) 


(77.4) 


sinr? 
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d*F = p. a {pr a + l)r^ a a a \ a sin ddr Add A dp + Va ^ Ua -Hi — — dt A dr A did , (£.5) 

sinri 

so the Yang-Mills equations (B. 3) (at J = 0) turn into the system 

Ma(Pa + l)r M “a a A a sin 2 d = -g 2 r fJ ’ a+Ub+Uc a a /3 b (3 c [[\ a , A ft ], A c ] , (E. 6 ) 

M"a - l)r"“- 2 /3“A a = — ff V-+^+Mc a a /3V[[A a , A 6 ], A c ] , (£.7) 

g{da - Pb)r^ a+ ^ b a a a b [\ a , A & ] = 0 , g(u a - v b )r Va+Vb (3 a (3 b [ A a , A 6 ] = 0 , a < 6, . (£.8) 


Then, obviously, in the left-hand sides of ( E.6-E.7 ) we should put all fi a = — 1 and all 
u 0 = 1 , respectively, for obtaining the confining solutions while the right-hand sides of 
(E.6-E.7) arc identically equal to zero only if \„j KC = A;j^. Under the circumstances 
the equations (E. 8 ) arc satisfied automatically and we come to the same conclusions about 
uniqueness of the confining solutions as in Subsection 8.1. 
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